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EXECUTIVE  SUMMARY 


Many  methods  have  been  proposed  for  suppressing  interference  in  a  direct  sequence  spread  spectrum 
signal  and  many  have  been  shown  to  be  effective  for  particular  jamming  conditions.  However,  it  has  been 
difficult  to  compare  the  performance  of  these  techniques  because,  in  most  cases,  they  have  been  studied 
separately  using  different  channel  conditions.  The  objective  of  this  effort  was  to  perform  a  direct 
comparison  of  a  number  of  techniques  by  evaluating  their  performance  for  identical  interferers.  The 
particular  techniques  under  consideration  were  predictive  and  two-sided  transversal  adaptive  filters,  the 
time  domain  adaptive  correlator,  transform  domain  excision  using  modulated  lapped  transforms  and 
transform  domain  detection,  and  transform  domain  excision  using  modulated  lapped  transforms  and  time 
domain  correlation.  The  interferers  included  fixed  and  variable  firequency  tones,  swept  tones  (chirps), 
narrowband  Gaussian  interference,  pulsed  tones  and  pulsed  narrowband  Gaussian  interference.  This 
collection  of  interferers  is  sufficiently  varied  to  test  the  ability  to  each  technique  to  resolve  interference  in 
both  fi-equency  and  time.  Additionally,  in  comparing  the  techniques,  the  computational  complexity  was 
taken  into  account,  such  that  direct  comparison  is  made  between  interference  suppression  techniques  of 
similar  complexity.  In  particular,  whenever  possible,  the  parameters  of  the  interference  suppressor,  such 
as  the  order  of  the  adaptive  filter  and  length  of  the  transform,  were  selected  to  make  the  complexity  nearly 
equivalent.  The  comparisons  are  based  on  bit-error-rate  performance  obtained  through  Monte  Carlo 
simulation.  Simulations  were  performed  using  the  Signal  Processing  WorkSystem  (SPW)  from  AltaGroup 
of  Cadence  Design  Systems,  Inc. 

The  simulation  results  indicate  that  the  performance,  in  terms  of  low  bit  error  rate,  of  the 
different  receiver  systems  depends  heavily  on  the  jamming  signal  and  parameters  used,  and  therefore 
there  is  no  “ideal”  structure  in  the  sense  of  a  particular  system  that  would  be  able  to  handle  any  type  of 
situation.  Nevertheless,  the  LMS  filter  systems,  'particularly  the  two-sided  filter,  present  excellent 
performance  in  removing  the  interference  when  the  channel  introduces  fixed-tone  jamming,  while  the 
tracking  capabilities  of  the  MLT  with  transform  domain  excision  and  detection  structure  are  remarkable. 
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CHAPTER  1 
Introduction 


Spread  spectrum  techniques  were  initially  conceived  for  military  usage,  due  to  their  ability  to 
reduce  the  effects  of  strong  jamming  and  provide  some  level  of  security.  However,  civilian  applications 
may  also  benefit  from  it,  such  as  the  use  of  this  system  for  mobile  radio  networks  (radio  telephony,  packet 
radio,  and  amateur  radio),  and  multiple-access  communication  (CDMA). 

The  processing  gain  of  a  direct  sequence  spread-spectrum  system  provides  some  degree  of 
protection  against  narrowband  interference.  However,  when  the  processing  gain  is  unable  to  combat  the 
jammer  adequately,  additional  signal  processing  techniques  can  be  employed  to  more  effectively  remove 
the  interference. 

Considerable  amount  of  work  has  been  performed  on  interference  mitigation  techniques 
involving  time  domain  and  frequency  domain  processing.  The  time  domain  study  has  been  centered 
primarily  on  transversal  filter  structures,  and  the  various  adaptive  algorithms  that  can  be  used  with  them. 
Transform  domain  processing  has  moved  from  initial  work  in  frequency  domain  excision  to  include  the 
consideration  of  different  transforms.  The  purpose  of  this  study  is  to  accomplish  a  comparative  study  of 
different  interference  suppression  techniques,  and  determine  which  technique  is  best  suited  for  a  given 
interference  scenario. 


1.1-  Spread  Spectrum  Systems 

Spread  spectrum  systems  transmit  the  data  using  a  much  larger  bandwidth  than  the  minimum 
bandwidth  required  [1, 2].  The  spectrum  spreading  is  achieved  before  transmission  employing  a  noise-like 
spreading  code  called  pseudo-random  or  pseudo-noise  sequence,  that  is  independent  of  the  message  sent. 
In  the  receiver,  the  same  code  is  used  to  despread,  or  correlate,  the  received  data,  which  allows  for  the 
recovery  of  the  original  message.  The  widening  of  the  bandwidth  is  realized  through  the  use  of 
modulation,  which  can  be  achieved  using  a  variety  of  methods,  including  pseudo-noise  (PN).  frequency 
hopping  (HP),  time  hopping  (TH),  and  hybrid  combinations  of  these  techniques.  In  pseudo-noise  direct- 
sequence  modulation,  a  pseudo-randomly  generated  sequence  is  used  to  produce  phase  transitions  in  the 
carrier  containing  the  data.  Frequency  hopping  involves  a  shift  in  the  carrier  s  frequency  in  a  pseudo¬ 
random  way,  while  in  time  hopping,  bursts  of  signal  are  originated  at  pseudo-random  times.  In  all  these 
techniques,  the  idea  behind  the  term  pseudo-random  is  random  in  appearance,  but  that  can  be  reproduced 
by  authorized  receivers. 
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An  important  attribute  of  a  spread-spectrum  communication  system  is  its  capability  to  reject 
externally  generated  interfering  (jamming)  signals  with  finite  power,  whether  intentional  or  not.  The 
transmitted  signal  assumes  a  wideband  noise-like  appearance,  and  therefore  is  hard  to  detect  from  the 
background  noise. 

A  pseudo-noise  (PN)  sequence  is  defined  as  a  random  binary  coded  sequence,  usually  periodic,  of 
Is  and  Os  with  certain  autocorrelation  properties.  Such  sequences  can  be  generated  by  simple  mechanisms, 
particularly  a  linear  feedback  shift  register.  Maximum-length  sequences  present  many  of  the  properties 

observed  in  truly  random  binary  sequences.  They  are  always  periodic  with  a  period  of  iV  =  2"  —  1 ,  where 
m  is  the  length  of  the  shift  register.  In  an  augmented  PN  sequence,  one  additional  0  bit  is  generated, 

thereby  increasing  the  period  to  N  =  2" .  Typically,  each  bit  interval  is  subdivided  into  N  short  intervals 
called  chips.  The  longer  the  period  of  the  spreading  code,  the  closer  will  the  transmitted  signal  be  to  a 
truly  random  wave,  and  the  harder  it  will  be  to  detect.  However,  the  price  for  a  better  protection  against 
interference  is  increased  transmission  bandwidth,  system  complexity,  and  processing  delay. 

In  (Mrder  to  better  illustrate  the  previous  concepts,  the  PN  direct  sequence  spread-spectrum 
modulation  operation  is  described  in  the  remainder  of  this  section: 

The  transmitted  message  b{t)  is  multiplied  by  a  wideband  PN  waveform  c{t) , 


m(t)  =  b(0  C(t).  (1.1) 

Thus,  if  the  data  is  narrowband  and  the  PN  sequence  is  wideband,  the  product  signal  /n(t)  will  have  its 
spectrum  bandwidth  about  the  same  as  the  PN  sequence,  as  it  equals  the  convolution  of  the  specu’a  of  the 
two  component  signals.  If  it  is  assumed  that 


(1.2) 


where  Tj,  is  the  bit  duration,  Ij.  is  the  chip  duration,  is  the  energy  per  symbol,  pit)  is  a  basic  pulse 

Tt 

shape,  and  b^  the  binary  digits.  iV  =  —  is  the  number  of  chips  per  data  bit. 

It  follows  that  the  power  spectral  density  of  the  message  bit)  is  given  by 


^  |'’(/)r  t, «»»(»)  .  (‘-3) 

n=— « 

and  assuming  that  its  autocorrelation  function,  R/,/,  in)  =  5, ,  the  delta  function, 

S,,if)  =  T,smc^itt/T,).  (1.4) 

and  that  c(f)  and  bit)  are  independent: 

S^i/)  =  SM)®S,,i/),  (1.5) 
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where  S^if)  and  are  the  power  spectral  densities  of  the  signals  tn{t)  and  c{t) 

respectively. 

For  N  »  I,  5„(/)  is  much  wider  than  therefore 

The  received  signal  is  modeled  as 

r(t)  =  m{t)  +  n(t)  + 1(0  (1-7) 

where  nif)  is  wideband  noise  and  i(f)  a  narrowband  interference  component.  To  recover  the  data  6(f) , 
the  received  signal  r(f)  is  applied  to  a  demodulator,  employing  the  same  PN  sequence  used  in  the 
transmitter  and,  assuming  that  the  receiver  operates  in  perfect  synchronism  with  the  transmitter, 

z(t)  =  c(t)  r(t) 

=  c^(t)  b(t)  +  c(t)  n(t)  +  i(t)  c(t) 

=  b(t)  +  n’(t)  +  i(t)  c(t),  (1.8) 

where  C^(t)  =  1  and  nXt)  is  a  new  wideband  noise  component.  If  n(f)is  zero-mean  WGN,  then 

n'{t)  is  likewise  a  zero-mean  WGN  process  with  double-sided  noise  spectral  density  ^  Watts/Hz. 

Similarly,  the  despread  interference  component  is 

i'(0  =  c(f)/(0 

=  (1.9) 

k 

where  represents  the  binary  digits. 

The  multiplication  operation  brings  the  data  signal  down  to  its  original  bandwidth  while  spreading  the 
noise  and  interference  to  a  wide  bandwidth,  as  shown  in  Figure  1.1. 
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Figure  1.1:  Slustration  of  Power  Spectral  Density  for  Spread  Spectrum  System 


Subsequent  lowpass  filtering  may  reject  large  amounts  of  the  spread  noise  and  interference, 
making  possible  the  recovery  of  the  original  message.  At  the  receiver,  the  ratio  of  the  spread  signal-to- 
noise  ratio,  SNRf ,  to  the  signal-to-noise  ratio  after  despreading,  SNR^ ,  is  defined  as  the  processing 
gain  Gp ,  i.e.. 


G 


p 


SNR, 
SNR,  ’ 


(1.10) 


and  is  a  measure  of  the  amount  of  performance  improvement  that  is  achieved  through  the  use  of  spread 
spectrum.  The  processing  gain  can  be  approximated  by  taking  the  ratio  of  the  spread  bandwidth  to  the 
data  information  bandwidth  (rate),  i.e.. 


G 


p 


(1.11) 


If  i(t)  is  orthogonal  to  c(t)  then  after  subsequent  integration  (matched  filtering),  the 
interference  component  can  be  eliminated.  Otherwise,  i'(t)  is  a  wideband  component  whose  effect  after 
subsequent  low-pass  filtering  is  the  same  as  if  i'(t)  were  Gaussian  with  power  spectral  density 

S.r(f)  =  S,(f)<8)SM) 
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N 


(1.12) 


Assuming  that  i(?)  is  narrowband  relative  to  c(t)  ,  with  Ef  the  interference  energy  in  a  bit  interval, 
and  and  S^{f)  the  power  spectral  densities  of  the  signals  and  i(t)  respectively. 

Finally 

n”(t)  =  n'(t)  +  i'(t)  (1.13) 

can  be  considered  Gaussian  with  power  spectral  density 

2E, 


2  2  N  2 


1  + 


NN„ 


(1.14) 


For  binary  antipodal  signaling,  the  bit  error  probability  is 


/ 

2E, 

2E,  ' 

1 

N, 

1  + - — 

111 

o 

NNg 

/ 

(1.15) 


where  is  the  energy  per  bit.  Suppose  that  the  jammer  has  average  energy  E,  but  only  on  a  fraction 
p  of  a  bit  interval.  The  noise  spectral  density  is  then 

N, 


(Vq  Ejlp 

2  N 


(1.16) 


when  the  jammer  is  on,  and 


No  N, 


.  2  2 

when  the  jammer  is  off.  The  error  probability  is  then 

f 


(1.17) 


r 


P,  =  il-P)Q 


2E, 


N 


+  pQ 


0  ) 


2E, 

N„ 

L  2£,/pl 

L  NNg  i 

(1.18) 


where  2  (jc)  = 


Suppose  that  Ej  /  N  »  Ng  /  2 ,  i.e. 
Then 


is  named  the  tail  function. 

the  processing  gain  is  insufficient  to  render  interference  negligible. 
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(1.19) 


P  Q 


pEy 
E,IN) 


An  intelligent  jammer  will  maximize  this  by  choice  of  p .  The  maximum  value  of  X  Q  {“Jx^  can  be 
seen  at  =  1.44  with  maximum  value  0.1657  .  The  maximum  value  of  p  is  then 


p*  =  min  j  1 ;  ■ 


1.44 


EylE;\ 

where  E'j  =  Ej  !  N  .  The  corresponding  error  probability  is 

,  0.1657 

‘  ~  Ey  IE',' 

which  decreases  only  algebraically  in  Ey  I E',. 

For  full-time  jammer,  the  error  probability  is 


/ 

2E, 

<s 

1 

1-1- — ^ 

ll 

L  ^o\ 

y 

and  if  E',  =  E ,  I N  »  N ^  1 2 ,  then 


Py-Q 


.  1 

<  —  e  ^ 

2 

which  decreases  exponentially  in  E,^  I E',. 


(1.20) 


(1.21) 


(1.22) 


(1.23) 

(1.24) 
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1.2-  Report  Overview 


An  elementary  way  to  fight  noise  and  interference  is  to  increase  the  signal  power.  However, 
raising  power  is  costly  and  can  damage  instruments.  Noise  modulation,  i.e.  utilization  of  noise-like  signal 
as  a  carrier,  is  very  suitable  for  anti-jamming  and  security  issues.  Nevertheless,  despite  the  fact  that 
spread-spectrum  modulation  techniques  inherently  present  interference  suppression  capabilities,  these  can 
be  enhanced  considerably  by  processing  the  received  signal. 

The  purpose  of  this  report  is  to  perform  a  direct  comparison  of  time  domain  and  firequency 
domain  interference  suppression  techniques.  Simulations  will  be  performed  using  the  Signal  Processing 
WorkSystem  (SPW)  from  AltaGroup  of  Cadence  Designs,  Incorporated.  The  overall  goal  will  be  to 
determine  the  performance  of  a  number  of  interference  suppression  techniques  under  different  interference 
conditions.  The  four  individual  techniques  studied  are  time  domain  adaptive  transversal  filtering, 
transform  domain  excision,  time  domain  adaptive  correlator  receiver,  and  transform  domain  excision  and 
detection.  A  direct  head-to-head  comparison  of  the  performance  of  the  techniques  for  a  set  of  “typical” 
interference  types  will  been  performed,  in  order  to  decide  which  technique  is  best  for  different 
applications.  The  types  of  interference  that  will  be  considered  are  stationary  tone  jamming,  swept  tone  and 
chirp  jamming,  narrowband  Gaussian  jamming,  pulsed  tone  and  pulsed  Gaussian  jamming.  The 
performance  measures  that  will  be  used  include  the  bit-error-rate  (BER),  tracking  performance,  and 
implementation  complexity. 

In  chapter  2,  the  different  systems  are  described,  beginning  with  the  modulator  and  channel, 
followed  by  the  various  jammer  types,  and  finally  the  four  communication  receiver  structures  are 
explained.  Besides,  the  implementation  complexity  for  the  receiver  systems  is  calculated.  In  chapter  3,  the 
Signal  Processing  WorkSystem  (SPW)  and  the  three  performance  measures  are  presented,  along  with  a 
description  of  the  simulations  to  be  performed.  In  chapter  4,  the  results  of  the  simulations  executed  for 
each  of  the  structures  with  the  different  types  of  jammers  are  presented,  along  with  the  curves  showing  the 
bit  error  rates  obtained.  Finally,  in  chapter  5,  the  discussions  and  conclusions  of  this  report  are  exposed, 
making  a  direct  comparison  among  the  different  interference  suppression  techniques  analyzed. 
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CHAPTER  2 
System  Description 


In  the  following  chapter,  the  different  systems  built  in  this  report  are  described. 

2.1-  Modulator  and  Channel 

A  block  diagram  showing  the  modulator  and  channel  (without  jamming)  can  be  seen  in  Figure 

2.1.  Binary  Phase  Shift  Keying  (BPSK)  baseband  signal  and  direct  sequence  spread-spectrum  modulation 
are  used  for  this  system. 

The  pseudo-noise  (PN)  sequence,  after  conversion  of  Os  and  Is  to  -Is  and  Is,  is  shown  in  Figure 

2.2.  A  value  of  6  is  used  for  the  length  of  the  shift  register,  with  a  resulting  period  of 

=  2"  —  1  =  2*  —  1  =  63 .  The  random  data  consists  of  random  Is  and  Os  binary  bits  (see  Figure  2.3), 
whose  sampling  frequency  equals  the  period  of  the  PN  sequence,  and  bit  rate  is  set  to  1.0  Hz .  The  logical 
Os  and  Is  are  converted  to  real-number  outputs  for  both  sequences:  when  the  input  is  true  (high;  >  0),  the 
output  is  set  to  1.0;  when  the  input  is  false  (low;  <  0),  the  output  is  set  to  -1.0.  The  PN  sequence  is  then 
used  to  modulate  the  random  data.  A  Gaussian  random  noise  (see  Figure  2.4)  is  added,  using  a  value  of 
Ef^  !  Nq  =  6.0  dB  for  the  following  parameters, 

Mean  =  0.0  (2.1) 

101ogio^Y^-60 

=Variance  =  10  =10  =7.91244.  (2.2) 

Table  2.1  summarizes  the  main  parameters  used  to  design  the  communication  channel. 


1  PN  Sequence 

m  =  Shift  Register  =  6  1 

N  =  Period  =  63  1 

Chip  Rate  =  63  Hz  j 

Random  Data 

Sampling  Frequency  =  63.0  Hz  | 
Bit  Rate  =  1.0  Hz  1 

Gaussian  Random  Noise 

Eh/N„  =  6.0  dB 

p  =  Mean  =  0.0 

=  Variance  =  7.91244 

Table  2.1:  Summary  of  Main  Parameters  for  Modulator  and  Channel 
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Dde  Order 


PN  SEQUENCE 


Figure  2.2:  PN  Sequence  Waveform 
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RANDOM  DATA 


Figure  2.3:  Random  Data  Waveform 
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NB  GAUSSIAN 


Figure  2.4:  Gaussian  Random  Noise  Waveform 
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2.2-  Jammers 


Intaference  is  unwanted  energy  from  systems  other  than  the  transmitter.  It  occurs  in  radio 
communication  when  the  receiver  gets  two  or  more  signals  in  the  same  frequency  band;  however,  it  may 
also  come  from  multipath  propagation  or  electromagnetic  transmission  cables.  Ctmtaminating  noise 
normally  has  an  additive  effect,  meaning  that  noise  adds  to  the  information-bearing  signal  at  different 
points  in  the  system.  The  sources  of  contamination  are  generally  classified  as  human  interference  or 
naturally  occurring  noise,  like  the  unavoidable  thermal  noise  [2]. 

In  this  report,  four  different  types  of  jammers  are  considered,  and  can  be  seen  in  figure  2.5 

[2,13]. 


2.2.1-  Stationary  Tone  Jamming  (Jammer  1).  A  sinusoid  of  fixed  frequency  is  used  as  the  jammer 
(Figure  2.6).  It  is  the  easiest  of  all  jamming  signal  to  generate,  and  can  be  modeled  as 

j[n]  =  Ay  cos[ 2 ;r  fjn  +  6  ],  (2.3) 

with  its  power  defined  relative  to  that  of  the  desired  signal  by  the  jammer-to-signal  ratio  (JSR) 


JSR  = 


(2.4) 


The  sampling  frequency  is  set  to  a  normalized  value  of  1.0  Hz ,  and  the  amplitude  (power)  of  the  jammer 


is  calculated  as 


JSR  M 

Aj  =  1.414  X 10  20  =  1.414  X 1020  =  14.14. 

Both  the  frequency  and  strength  of  this  jammer  are  parameters  to  be  varied. 


(2.5) 


2.2.2-  Swept  Tone  and  Chirp  Jamming  (Jammer  2).  In  the  case  of  a  swept  tone,  the  frequency  of  the 
sinusoid  is  varied  using  a  triangular  wave  (Figure  2.7),  while  for  chirp  jamming  the  frequency  is  varied 
using  a  ramp  wave  (Figure  2.8).  Here,  the  jammer  sweeps  through  the  frequency  band  of  the  transmitted 
signal,  meaning  that  its  frequency  is  now  constandy  changing,  and  therefore  the  communication  receiver 
systems  must  also  be  constantly  changing  to  track  the  jammer.  This  jammer  can  also  be  represented  by 
equation  (2.3),  with  a  moving  value  for  fj .  The  value  of  JSR  can  also  be  calculated  according  to 
equation  (2.4),  while  the  sweep  frequency  should  allow  the  jammer  frequency  to  sweep  back  and  forth 
many  times.  The  frequency  range,  variation  rate  and  power  are  the  parameters  to  be  varied  for  this  type  of 
jammer. 
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JAMMER  1  -  Fixed  tone 


Figure  2.6:  Stationary  Tone  Jamming  Waveform 
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JAMMER  2  TRIANG  -  Swept  Tone 


Figure  2.7:  Swept  Tone  Waveform 
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Figure  2.8:  Chirp  Jamming  Waveform 


3000 
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2.2.3-  Narrowband  Gaussian  Jamming  (Jammer  3).  The  output  from  a  white  noise  generator  with  a 
Gaussian  distribution  is  passed  through  a  narrowband  filter  to  produce  the  interference  (Figure  2.9).  This 
class  of  jammer  is  a  brute-force  jammer  that  does  not  exploit  any  knowledge  of  the  antijam 
communication  system  except  for  its  spread  bandwidth. 

Denoting  the  total  interference  power  as  ,  the  jammer-to-signal  ratio  is  given  by 


JSR  = 


^nb 


(2.6) 


The  filter  is  an  HR  Butterworth  bandpass  filter,  of  order  10  and  attenuation  at  passband  edge  3  dB , 
while  the  sampling  frequency  is  set  to  a  normalized  value  of  1.0  Hz .  The  main  parameters  for  this 
jammer  are  calculated  as  follows, 

/I  =  Mean  =  0.0  (2.7) 

JSK 

2  10 

a  =  Variance  = - .  (2.8) 

2  BWj  ^  ’ 

Filter  bandwidth  and  jammer  power  are  the  parameters  to  be  varied. 


2.2.4-  Pulsed  Tone  and  Pulsed  Gaussian  Jamming  (Jammer  4).  A  pulse-shaped  masking  function  is 

placed  on  tone  (Figure  2.10),  and  Gaussian  (Figure  2.11)  jammers.  The  sampling  frequency  is  set  to  a 

normalized  value  of  1.0  Hz ,  and  the  tone  amplitude  Aj  =  14.14  is  calculated  according  to  equation 

(2.5).  As  to  the  Gaussian  noise,  the  filter  is  an  HR  Butterworth  bandpass  filter,  of  order  10  and 

attenuation  at  passband  edge  3  dB  .  The  masking  function  consists  of  a  square  wave  varying  from  0  to 
■> 

1 .  The  value  of  the  quiescent  frequency,  i.e.,  the  frequency  of  the  masking  function,  has  to  be  lower  than 
1/ Period =  l/A^  =  l/63,  and  is  set  to  0.000001  Hz.  The  duty  cycle  of  the  masking 

function  is  a  parameter  that  specifies  the  percentage  of  time  that  the  signal  is  high,  rather  than  low,  and  is 
set  to  Factor^^  =  50  %  . 

The  main  parameters  for  this  jammer  are  calculated  as  follows. 


fi  =  Mean  =  0.0 

(2.9) 

JSR 

10 

<7^  =  Variance  = - . 

2 BW,  Factor.^ 

(2.10) 

The  parameters  of  the  original  tone  and  Gaussian  jammers  are  quantities  to  be  varied. 
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JAMMER  4  -  TONE  w/  Mask  Func 


Figure  2.10:  Pulsed  Tone  Jamming  Waveform 


21 


2.2.5-  Summary 

Table  2.2  gives  the  default  parameter  values  used  for  the  different  jammers.  Unless  otherwise 
indicated,  these  are  the  values  that  are  used. 


JAMMER  1 

Eb/No  =  6dB  I 

JSR  =  20.0dB 

Jammer  Frequency  =  0.234  Hz 

JAMMER  2 

Eb/No  =  6dB 

JSR  =  20.0dB 

Sweep  Frequency  =  0.00001  Hz 

JAMMER  3 

Eb/No  =  6dB 

JSR  =  20.0dB 

Jammer  Frequency  =  0.234  Hz 

Jammer  Bandwidth  =  0.1  Hz 

JAMMER  4 

Eb/No  =  6dB 

JSR  =  20.0dB 

Jammer  Frequency  =  0.1  Hz 

Jammer  Bandwidth  =  0.05  Hz 

Masking  Function  Frequency  =  0.000001  Hz 

Masking  Function  Duty  Cycle  =  50% 

Table  2.2:  Summary  of  Main  Parameters  for  Jammers 


2.3-  Communication  Receivers 

The  communication  receivers  detect,  demodulate  and  despread  the  received  data  signal.  Besides, 
they  perform  additional  processing  in  an  effort  to  reconstruct  or  estimate  the  originally  transmitted 
message  symbol,  contaminated  by  the  jamming  process. 

2.3.1-  Adaptive  Filters 

Filters  are  used  in  nearly  any  communication  system,  for  the  purpose  of  retrieving  an 
information-bearing  signal  from  unwanted  contaminations  such  as  interference,  noise,  and  distortion 
products.  When  the  correlation  functions  are  unknown  and  the  environment  is  non-stationary,  adaptive 
filters  can  be  an  excellent  alternative.  These  filters  are  self-designing  devices  whose  coefficients  are 
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updated  according  to  a  recursive  algorithm,  such  that  efficient  transmission  is  maintained.  The  LMS 
(Least-Mean-Squared)  algorithm  is  probably  the  best  known  of  the  adaptive  filtering  algorithms.  A 
significant  feature  of  the  LMS  algorithm  is  its  simplicity.  Indeed,  this  algorithm  is  iterative,  i.e.  it  does  not 
require  information  about  the  statistics  of  the  incoming  signal,  but  rather  adapts  on  a  step-by-step  basis  to 
converge  to  the  optimal  (in  the  mean-square  sense)  tap  weights.  In  the  direct  sequence  receiver 
application,  an  iteration  occurs  each  time  a  new  sample  is  input  into  the  delay  line,  meaning  that  the 
system  adapts  once  per  chip  of  the  signal  [3, 4, 5]. 

Since  both  the  data  and  the  thermal  noise  are  wideband  signals,  their  future  values  cannot  be 
easily  predicted  from  their  past  values.  However,  the  interference  is  a  narrowband  process,  and  hence  can 
have  its  future  values  predicted  from  past  values.  Once  its  present  value  is  predicted,  it  can  be  subtracted 
from  the  received  signal.  The  same  principle  is  used  for  the  two-sided  transversal  filter,  except  that  now 
the  estimate  of  the  current  value  of  the  interference  is  based  upon  both  past  and  future  values. 

In  figure  2.12  ,  a  simple  predictive  adaptive  filter  interference  suppression  structure  is  shown. 
The  received  signal  x[k]  is  fed  into  a  tapped  delay  structure  having  a  tap  spacing  equal  to  the  chip 
duration,  .  Delayed  input  samples,  —  1] ,  X^k  —  2] , x[k  —  A^]  are  multiplied  by  a  set  of  tap 
weights,  Wp  Wj  >  •••»  .  and  the  products  are  summed  to  form  the  signal  y[^].  This  output  is  then 

subtracted  from  the  input  signal  to  form  an  error  signal,  €  [^] .  The  tap  weights  are  determined  using  the 
LMS  adaptive  algorithm  which  works  to  minimize  the  error  signal,  e  [^] .  This  adaptive  process  forces 
y{,k\  to  become  an  estimate  of  the  current  sample  of  the  input  signal,  x{k^ . 

The  LMS  algorithm  for  adaptive  filters  is  as  follows 

w-Xk  +  l]  =  wXk]+^e  [A:]  x[k-i] .  (2.11) 

The  term  /i  e  [it]  x\^k  —  i]  represents  the  correction  that  is  applied  to  the  current  estimate  of  the  tap- 
weight  vector,  .The  adaptation  constant  /i  is  a  parameter  which  determines  the  rate  of  convergence 
of  the  algorithm.  Indeed,  the  LMS  algorithm  is  convergent  in  the  mean  square  if 

0<u  < - - - .  (2-12) 

Total  Average  Input  Power 

The  larger  the  adaptation  constant  p  ,  the  faster  the  tracking  capability  of  the  LMS  algorithm,  which  is 
equivalent  to  the  LMS  algorithm  having  a  short  “memory”.  However,  a  large  value  of  p  may  imply  a 
high  excess  mean-square  error,  sometimes  called  misadjustment,  as  limited  data  is  applied  to  estimate  the 
gradient  vector.  Therefore,  a  compromise  must  be  reached  between  fast  tracking  and  low  excess  mean- 
squared  error.  Hence,  the  parameter  p  may  be  conceived  as  the  memory  of  the  LMS  algorithm,  since  it 

decides  on  the  weighting  applied  to  the  tap  inputs. 
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Figure  2.12:  Predictive  Transversal  Filter 
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The  inputs  to  the  predictive  filter  are  a  delayed  version  of  the  signal,  coming  from  the  delay 
element,  and  the  desired  response,  which  is  the  signal  coming  directly  from  the  channel.  The  filter  length 
for  the  predictive  LMS  filter  system  is  set  to  7  taps.  In  the  two-sided  system,  the  signal  is  fed  into  a  first 
LMS  adaptive  filter.  The  output  of  the  second  filter,  whose  inputs  are  a  delayed  version  of  the  signal  and 
the  error  feedback,  is  summed  up  with  the  output  of  the  first  filter  and  subtracted  from  the  output  of  the 
first  unit  delay  to  fcffm  the  error  signal  for  both  LMS  filters.  The  filter  length  for  the  two-sided  LMS 
adaptive  filter  system  is  set  to  5  taps. 


2.3.2-  Modulated  Lapped  Transforms  -  Transform  Domain  Excision  (MLT-TDE) 

The  basic  idea  behind  transform  domain  processing  is  to  compactly  represent  the  interference 
energy  in  the  transform  domain,  and  then  easily  remove  it.  The  transform  should  be  selected  such  that  the 
interference  is  nearly  a  delta  function  in  the  transform  domain,  while  the  desired  signal  is  converted  to  a 
waveform  whose  spectrum  is  notably  wide  in  comparison  to  the  transformed  interference  s.  A  simple 
exciser  can  then  remove  the  interferer,  without  significantly  damaging  the  original  signal.  An  inverse 
transform  then  recovers  the  almost  interference-free  desired  signal  [6,  7,  8].  Figure  2.7  illustrates  the 
excision  process. 

The  excision  process  is  a  fast  adaptive  algorithm  that  is  performed  prior  to  the  correlation  stage. 
It  removes  narrowband  interference,  adjusting  the  position  of  the  notch  (or  notches)  to  suppress  spectral 
components  with  large  amount  of  energy,  typically  those  that  surpass  a  user-defined  threshold.  Its  level 
can  be  established  sufficiently  high  so  that  components  are  removed  only  when  interference  is  present. 
However,  if  the  interference  is  changing  rapidly,  the  algorithm  may  not  be  able  to  track  these  variations. 
The  performance  of  the  excision  process  highly  depends  on  the  ability  of  the  transform  to  compactly 
represent  the  interference. 
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(a)  Signal  +  Interference 


(c)  Excised  Transform 


Figure  2.13:  DIustration  of  the  Excision  Process 
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The  Fast  Fourier  Transform  (FFT)  is  not  considered  here,  due  to  its  considerable  spectral  overlap 
effects  with  large  number  of  adjacent  carriers.  Instead,  Lapped  Transforms  (LTs)  and  in  particular 
Modulated  Lapped  Transforms  (MLTs),  are  proposed  [7].  The  central  reason  for  the  introduction  of  LTs  is 
found  in  one  of  the  major  disadvantages  of  traditional  block  transforms:  the  blocking  effects,  which  are 
discontinuities  in  the  reconstructed  signal.  Besides,  LTs  have  higher  coding  gains  than  block  transforms, 
and  therefore  higher  signal-to-noise  ratios.  The  basis  functions  of  LTs  are  filters  whose  length  L  is  equal 
to  some  even  integer  multiple  of  the  number  of  filter  bank  subbands,  i.e.  Li  'IKM  ,  where  Af  is  the 
number  of  subbands  and  K  is  the  overlapping  factor.  MLTs  are  considered  as  a  subset  of  general  LTs 
with  K  =  \,  and  they  closely  approximate  the  ideal  filter  system  of  unity  passband  gain  and  infinite 
stopband  attenuation.  A  2M-tap  lowpass  filter  is  used  as  a  subband  filter  prototype  and  shifted  in 
frequency  to  produce  a  set  of  OTthogonal  bandpass  filters,  which  span  the  frequency  domain.  The  lowpass 
prototype,  ,  must  satisfy  the  following 

hi2M-l-n\  =  h[ri\  0<n<M-l  (2.13) 


and 


h\n]  +  h\n+M]  =  \  0<m<M/2-1. 

The  half-sine  windowing  function  is  used  as  the  lowpass  filter  prototype. 


/z[«]  =  -  sin 
The  basis  functions  can  then  be  described  as 


n  +  - 

\  2 


It 

2M 


Yk\-n]  =  h[n] 


[T 

,  —  cos  n+- 

VM  I 


M  +  lV.  1 


7t 


2)  M 


where  0 < /z ^ 2M  —  1  and  — 


(2.14) 


(2.15) 


(2.16) 


In  the  MLT-TDE  system,  the  narrowband  interference  is  removed  in  the  transform  domain.  The 
signal  is  brought  to  the  transform  domain  via  a  modulated  lapped  transform.  The  transform  domain 
excision  is  then  performed,  regulated  by  two  parameters,  the  threshold  and  the  notch  width.  The  exciser 
block  cancels  any  components,  and  a  user-defined  number  of  bins  around  the  components,  that  exceed  a 
preset  threshold.  Next,  the  signal  is  brought  back  to  the  time  domain  via  an  inverse  modulated  lapped 
transform.  The  vector  length  for  the  system  is  set  to  32  points. 

In  both  the  adaptive  filters  and  MLT-TDE  systems,  the  receiver  does  not  need  to  be  synchronized 
with  the  transmitter.  Both  receivers  can  be  seen  in  Figure  2.14. 
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2.3.3-  Adaptive  Correlator 


An  adaptive  correlator  is  a  system  that  performs  the  following  functions:  it  detects  the 
transmitted  data,  removes  interference,  and  compensates  for  multipath.  Consequently,  the  resulting 
receiver  would  consist  of  a  3  tapped-delay  line  structure  in  cascade:  the  LMS  adaptive  filter  to  remove 
interference,  the  correlator  to  despread  the  signal,  and  a  linear  combiner;  still,  a  single  fdter  may  be  used 
to  perform  all  three  tasks  [9, 10]. 

There  are  two  stages  of  operation: 

During  the  training  period,  a  known  data  sequence  is  transmitted  and  a  synchronized  version  of 
this  signal  is  generated  in  the  receiver  where  it  is  used  as  the  desired  response. 

When  the  training  is  completed,  the  adaptive  equalizer  is  switched  to  its  second  mode  of 
operation,  the  decision-directed  mode.  The  system  is  operated  by  aligning  a  data  bit  in  the  adaptive 
correlator,  summing  the  products  of  the  tap  weights  and  signal  samples,  and  comparing  the  result  to  the 
original  transmitted  data  bit.  The  error  between  the  output  of  the  adaptive  correlator  and  the  transmitted 
data  bit  is  then  used  to  adapt  the  tap  weights  according  to  the  LMS  algorithm,  in  order  to  minimize  the 
mean-square  error.  If  the  length  of  the  correlator  is  the  same  as  that  of  the  data  bit,  timing  has  to  be 
synchronized  between  the  receiver  and  the  transmitter,  and  the  adaptation  has  to  be  exactly  one  correlator 
tap  per  received  signal  sample,  since  the  adaptive  algorithm  is  accomplished  once  per  data  bit.  In  a  direct 
sequence  spread  spectrum  signal  in  AWGN,  the  adaptive  correlator  will  converge  to  the  standard 
correlator  receiver  by  in  fact  learning  the  transmitter’s  spreading  sequence  and  the  channel.  When 
interference  is  present,  the  adaptive  correlator  will  tend  to  suppress  the  interference.  However,  the 
adaptive  correlator  can  jump  to  the  wrong  polarity,  producing  high  error  rates  by  inverting  all  the  bits 
after  the  jump  occurs. 

Figure  2.8  is  a  block  diagram  of  the  adaptive  correlator  receiver.  The  receiver  is  a  fractionally- 
spaced  linear  equalizer  having  a  tap  spacing,  ,  that  is  some  fraction  of  the  chip  period  and  a  total 
duration,  ,  that  equals  the  data  bit  period  or  more.  During  the  operation  of  the  receiver,  the 

contents  of  the  delay  line  will  be  advanced  by  a  full  data  bit  period  between  iterations  of  the  adaptive 
algorithm.  As  a  result,  a  particular  tap  of  the  tapped  delay  line  will  always  correspond  to  a  received  signal 
sample  taken  at  the  same  time  within  the  data  bit  period. 
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Figure  2.15:  Adaptive  Correlator  Receiver  Structure 
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The  signal  is  fed  into  a  Least  Mean-Squared  (LMS)  adaptive  filter.  The  output  of  the  filter  is  then 
subtracted  from  the  original  signal  data  to  form  the  error  feedback  of  the  filter.  A  clock  holds  the  filter 
operation  for  a  period  of  63  iterations  (period  of  the  PN  sequence),  making  the  correlator  adapt  once  every 
data  bit.  The  filter  length  is  therefore  set  to  the  same  value  as  the  period  of  the  PN  sequence,  63  taps. 

The  signal  samples  in  the  delay  line  can  be  expressed  as  the  vector  y  [«] , 

y[«]=(  3'[«7;],}’[n7;  +(4  -i)rj)  (2.17) 

where  4  is  the  data  bit  period,  is  the  chip  period,  and  4  ‘S  the  length  of  the  spreading  sequence  in 
chips.  Likewise,  the  tap  weights  can  be  represented  as  a  vector. 

The  ouqjut  of  the  filter  d'{n\  is  found  by  taking  the  dot  product 

d\n\  =  y^[n\  vv*  =  vv^  t[/j]  .  (2.19) 

The  output  is  compared  to  the  desired  input  to  the  filter,  rf[n] ,  which  is  the  transmitted  data  bits  obtained 
during  the  training  period.  The  difference  between  d\n\  and  d[n\  is  the  error  signal  S ;  the  tap 
weights  of  the  filter  are  selected  to  minimize  its  value.  For  the  theoretical  analysis,  we  will  consider  a 
channel  containing  only  white  Gaussian  noise  and  a  single  tone  interferer.  The  power  at  the  receiver  due 
to  the  transmitted  signal  is  equal  to  P, ,  and  the  spreading  code  is  C .  The  received  vector  y{n\  is 


y[«]  =  X  [n]^  c  +nn]  +  n[n] .  (2.20) 

According  to  the  analysis  in  [10],  the  data  bit  energy  to  one-sided  noise  power  spectral  density  {Ej,  I N^) 


IS 


E,  _  4  Ps 


the  jammer-to-signal  power  ratio  (JSR)  is 


JSR  =  ^, 


and  the  minimum  mean  square  error  is 


e  =  ■ 


1  + 


PsA] 


(t"(ct"  +  4  A])\ 


(2.21) 


(2.22) 


(2.23) 


r  H 
c  e  e  c 


where  e  =  C  =  L. 
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2.3.4-  Modulated  Lapped  Transforms  -  Transform  Domain  Excision  and  Detection 
(MLT-TDED) 

In  the  MLT-TDED  system,  an  augmented  PN  sequence  is  used,  generating  a  period  of 
^  =  2"  =  64 ;  the  randcan  data  sampling  frequency  also  results  in  a  value  64.0  Hz .  All  the  vector 

lengths  are  set  to  64  points  in  this  system,  instead  of  the  63-length  used  in  the  three  previous  receivers. 

The  detection  is  performed  in  the  transform  domain  without  performing  an  inverse  transform.  Figure  2.9 
illustrates  a  generic  communications  receiver  employing  MLT-TDED.  In  this  system,  the  spectral 
coefBcients  of  the  observed  data  signal  are  point-wise  multiplied  with  the  weighting  vector  elements.  The 
ccHTesponding  outputs  are  then  cmrelated  with  the  known  transform  domain  reference  vectors,  using  a  dot 
product,  to  yield  the  decision  variable.  The  excision  process  is  performed  using  a  64  point-length  vector. 


Fot  both  the  adaptive  correlatw  and  MLT-TDED  receiver  systems,  timing  has  to  be 
synchronized  between  the  receiver  and  the  transmitter,  since  the  processes  are  accomplished  once  per  data 
bit.  Figure  2.17  presents  both  systems. 
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Figure  2.16:  Receiver  Employing  MLT  Domain  Excision  and  Detection 
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2.3.5-  Computational  Complexity 

In  Table  2.3,  the  computational  complexity  in  terms  of  OpS  (Operations  per  Sample)  is 

calculated  for  the  four  different  communication  receiver  systems. 

The  adaptive  filters  and  MLT-based  transform  exciser  are  configured  to  have  roughly  the  same 
complexity  in  order  to  make  a  fair  comparison  between  the  two  structures.  Therefore,  the  taps  of  the  filters 
and  the  vector  length  used  fw  the  MLT-TDE  system  are  selected  to  satisfy  this  requirement,  attaining 
complexities  of  28  OpS  fw  the  predictive  fitter,  25  OpS  for  the  two-sided,  and  29  OpS  for  the 

MLT-TDE  structure. 

On  another  hand,  the  adaptive  correlator  and  the  MLT  with  transform  domain  excision  and 
detection  are  constrained  by  the  need  to  have  the  block  size  (vector  length)  equal  the  length  of  the  data,  as 
the  receiver  requires  to  be  synchronized  with  the  transmitter  for  these  two  structures.  Values  of  4  OpS 
for  the  correlator  and  17  OpS  for  the  MLT-TDED  receiver  are  obtained. 


Predictive  LMS  Filter 

Two-Sided  Adaptive  LMS  Filter 

NpF  =  Filter  Order  =  7  taps 

OpS  =  4NpF  =  28 

Nuf  =  Filter  Order  =5  taps 

OpS  =  5NuF  =  25 

MLT-TDE 

Nmlt-tde  =  Vector  Length  =  32  points 

OpS  =  9  +  4  log2  Nmlt-tde  =  29 

j  Adaptive  CorrelatCM- 

Nc  =  Filter  Length  =  63  taps 

1 

OpS  =  2  +  2  =  4 

MLT-TDED 

Nmlt-tded  =  Vector  Length  =  64  points 

1 

OpS  =  5  +  2  log2  Nmlt-tded  =17 

Table  2,3:  Computational  Complexity  for  the  Different  Communication  Receiver  Systems 
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CHAPTER  3 
Implementation 


3.1-  Signal  Processing  WorkSystem  (SPW) 

The  Signal  Processing  WwkSystem  (SPW)  from  the  Alta  Group  of  Cadence  Design  Systems,  Inc. 
offers  the  possibility  of  capturing,  simulating,  testing  and  implementing  a  broad  range  of  digital  signal¬ 
processing  (DSP)  designs,  like  digital  communication  and  radar  systems.  Its  use  goes  from  the  evaluation 
of  various  architectural  approaches  to  the  design,  development,  simulation  and  fine-tuning  of  DSP 
algcvithms.  SPW  coisists  of  several  modules,  and  contains  different  libraries.  Besides,  it  can  access 
directly  DSP  boards  residing  in  the  computers  through  the  Code  Generation  System  (CGS),  which 
generates  C  code.  These  programs  are  then  automatically  downloaded  to  the  computer,  along  with  any 
required  files,  where  they  are  compiled  and  downloaded  into  the  DSP  board  for  execution.  This  process 
shortens  the  time  required  to  run  the  simulations. 


3.2-  Parameters 

The  three  primary  performance  measures  used  are: 

Bit-Error-Rate  (BER).  While  signal-to-noise  (SNR)  improvement  is  often  used  as  a  measure  of 
the  performance  of  interference  suppression  techniques,  BER  might  be  a  better  parameter,  since  it  takes 
into  account  both  the  level  of  interference  suppression  and  the  amount  of  distortion  introduced  into  the 
spread  spectrum  signal.  Fot  its  computaticxi,  the  received  data  signal  is  compared  with  the  original 
transmitted  data  signal,  from  which  a  bit  error  rate  is  generated,  according  to  the  following  equation. 


BER  = 


#  Erroneous  Bits 
#  Bits  Sera 


(3.1) 


The  theoretical  perfamance  fw  BPSK  in  AWGN  is  related  to  E^  I [dB]  as 

BER  =  erfc  Qe^  /  )  (3.2) 

where  erfc  is  the  complementary  enor  function  defined  as 

a» 

erfc(x)  =  ^^je-‘"dt.  (3.3) 

Implementation  Complexity.  The  implementatimi  cost  of  a  particular  technique  can  greatly 
detomine  its  applicability.  Digital  implementations  of  the  various  techniques  are  proposed,  and  their 
design  complexity  is  determined  in  terms  of  operations/chip  or  operations/bit. 

Tracking  Performance.  Many  interference  suppression  techniques  require  time  to  learn  the 
charact^stics  of  the  intaf(^ence  in  order  to  remove  it.  In  iterative  adaptive  systems,  this  learning  period 
is  often  called  Ae  convo-gence  time. 
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Clearly,  there  is  significant  interaction  between  the  various  performance  measures.  For  instance, 
the  BER  and  tracking  perframance  of  the  time  domain  adaptive  filter  is  a  function  of  the  number  of  taps 
which,  in  turn,  affects  the  implementation  complexity.  Performance  results  will  be  organized  in  a  manner 
which  allows  a  fair  comparison  between  the  various  techniques. 

3.3-  Simulations 

The  different  simulations  performed  for  the  four  jammer  systems,  using  the  previously  described 
communication  receivers,  can  be  seen  in  Table  3.1. 


Jammer  1 

BER  vs.  Eb/No 

BER  vs.  JSR* 

BER  vs.  Jammer  Frequency 

Jammer  2 

BER  vs.  Eh/No 

BER  vs.  JSR* 

BER  vs.  Sweep  Frequency 

Jammer  3 

BER  vs.  Eb/No 

BER  vs.  JSR* 

BER  vs.  Jammer  Frequency 

BER  vs.  Jammer  Bandwidth 

Jammer  4 

BER  vs.  Eh/No 

BER  vs.  JSR* 

BER  vs.  Jammer  Frequency 

BER  vs.  Jammer  Bandwidth 

Table  2 

(.1:  Simulations  Performed 

The  simulations  are  run  until  the  number  of  errors  obtained  reaches  the  value  of  200.0 .  There  is 
also  a  value  fw  the  number  of  symbols  to  wait  before  the  system  actually  starts  computing  the  BER.  This 
parameter  allows  the  system  to  stabilize  before  the  computation  of  the  BER  begins.  Its  value  is  set  to 
2,000  fOT  the  adaptive  transversal  filters,  MLT-TDE  and  MLT-TDED  systems,  while  a  value  of  20,000 
is  used  fw  the  adaptive  correlator  receiver. 


^  The  decision  parameters  for  the  different  receiver  systems,  namely  fJ.  for  the  LMS  algorithm, 
and  notch  and  threshold  for  the  excision  process,  are  optimized  for  each  value  of  JSR. 
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CHAPTER  4 
Results 

In  the  following  chapter,  the  curves  representing  the  Bit  Error  Rates  associated  with  the  four 
jammo's,  for  the  different  parameters  menticmed  in  Table  3.1,  are  shown. 

The  plots  are  presented  in  such  a  way  that  the  LMS  adaptive  filters  and  the  MLT-based  transform 
domain  excision  are  presented  in  the  same  figure,  while  the  LMS  adaptive  correlator  and  the  MLT-based 
transfcHm  domain  excision  and  detection,  are  displayed  together  in  another  figure.  This  division 
filters/MLE-TDE  and  corelatOT/MLT-TDED  is  done  in  order  to  perform  a  better  comparison  between 
systems  that  present  similar  characteristics.  Indeed,  the  LMS  adaptive  filters  and  the  MLT-TDE  system 
face  similar  implementation  complexity,  with  a  value  of  28  OpS  for  the  predictive  filter,  25  OpS  for 
the  two-sided,  and  29  OpS  fiw  the  MLT-TDE  structure  (refer  to  Table  2.3).  On  another  hand,  the 
ccOTelatcH'  and  MLT-TDED  systems  require  less  computation,  with  a  value  of  4  OpS  for  the  correlator, 
and  17  OpS  for  the  MLT-TDED  receiva.  Besides,  in  both  the  LMS  adaptive  filter  and  the  MLT-TDE 
systems,  the  receiver  does  not  need  to  be  synchronized  with  the  transmitter,  while  for  the  correlator  and 
MLT-TDED  structures,  timing  has  to  be  synchronized. 

4.1-  Results  for  Jammer  1 

This  section  presents  the  results  obtained  tor  each  of  the  receiver  structures  when  the  channel 
introduces  fixed-frequency  tone  jamming. 

Figure  4.1  plots  BER  vs.  Eh/No  fw  the  predictive  and  two-sided  LMS  filters  and  MLT-based 
transform  domain  exciser.  Also,  the  perfcmnance  of  the  BPSK  baseband  signal  in  AWGN  alone  (with  no 
jamming),  as  well  as  the  performance  obtained  without  any  interference  suppression  technique  (raw  error) 
are  shown.  Fot  these  simulaticxis,  the  default  parameters  used  for  the  jammer  are  JSR  =  20  dB  and  a 
normalized  frequency  f  =  0.234  Hz .  A  value  of  pL  —  0.0000001  for  the  predictive  and 
jl  =  0.00001  fw  the  two-sided  filtCT  are  used,  while  a  threshold  =  15.0  and  notch  width  =  1  are 
employed  fw  the  MLT-transfwm  domain  exciser.  These  parameters  were  optimized  in  order  to  obtain  the 
best  performance  for  the  receiver  systems.  As  can  be  seen  in  the  figure,  all  three  suppression  structures 
present  similar  p^cnmance,  and  a  significant  improvement  over  the  no-suppression  case. 
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Figure  4.2  shows  HER  vs.  Eh/N,  far  the  adaptive  LMS  correlator  and  MLT  with  transform 
domain  excision  and  detection.  Furthermore,  the  performance  of  BPSK  baseband  signal  in  AWGN  alone 
and  the  (me  obtained  without  any  interference  suppression  technique  (raw  error)  are  added.  For  these 
simulations,  the  default  parameters  used  for  the  jammer  are  JSR  =  20  dB  and  a  normalized  frequency 
/  =  0234  Hz  .  A  value  of  =  0.000(X)1  for  the  correlator,  and  threshold  =  25.0  and 
notch  width  =  1  the  MLT-TDED  system  are  used.  These  parameters  were  optimized  in  order  to 
achieve  best  performance  for  the  receivers.  Note  that  both  structures  present  similar  results,  and  that  their 
performance  is  closer  to  the  AWGN-only  results  than  the  previous  systems  presented  in  Figure  4.1. 
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JAMMER  1  -  Fixed  Tone 


Figure  4.2 ;  BER  vs.  Eb/No  for  Jammer  1  •  Corr  /  MLT-TDED 


Figures  4.3  and  4.4  show  the  curves  of  BER  vs.  JSR  for  the  different  structures,  using  a  value  of 
=  6  dB  fw  the  Gaussian  noise  in  the  channel  and  a  normalized  frequency  f  =  0.234  Hz  for 
the  jammer.  In  wdo"  to  obtain  these  results,  the  parameters  of  the  interference  suppression  techniques,  i.e. 
the  convergence  factw  fe  the  adaptive  filters  and  correlator,  and  exciser  parameters  for  the  MLT 

systems  are  optimized  for  each  value  of  JSR,  and  can  be  seen  in  Tables  4.1  and  4.2.  The  figures  indicate 
that  the  LMS  filters  provide  the  best  performance,  showing  little  dependence  of  BER  on  JSR  ,  obtaining 
better  results  with  the  two-sided  filter.  The  performance  of  all  the  other  systems  is  relatively  poor  after 
values  of  JSR  ~  30  dB,  though  the  performance  of  the  MLT-TDED  structure  is  slightly  better  than  the 

others. 
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Figure  4J  :  BER  vs.  JSR  for  Jammer  1  -  Filters  /  MLT-TDE 


Figure  4.4  :  BER  vs.  JSR  for  Jammer  1  -  Corr  /  MLT-TDED 
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Table  4.1:  Jammer  1  -  Filters  /  MLT-TDE  -  Parameters  Used  for  Each  Value  of  JSR 


ADAPTIVE  CORRELATOR 


MLT-TDED 


JSR=  10  dB 

JSR=10dB  1 

M.  =  0.000001 

Threshold  =  11.0 

JSR  =  15  dB 

Notch  Width  =  1 

^  =  0.000001 

JSR=  15  dB 

JSR  =  20  dB 

Threshold  =  11.0 

=  0.000001 

Notch  Width  =  1 

JSR  =  25dB 

JSR  =  20dB 

|i  =  0.000001 

Threshold  =  25.0 

JSR  =  30dB 

Notch  Width  =  1 

^1  =  0.000001 

JSR  =  25dB 

JSR  =  35dB 

Threshold  =  11.0 

H  =  0.000001 

Notch  Width  =  1 

JSR  =  40  dB 

JSR  =  30dB 

H  =  0.000001 

Threshold  =  13.0 

JSR  =  45dB 

Notch  Width  =  1 

|a  =  0.0000001 

JSR  =  35dB 

JSR  =  50  dB 

Threshold  =13.0 

^1  =  0.00000001 

Notch  Width  =  1  | 

JSR  =  55dB 

JSR  =  40  dB 

Threshold  =  13.0 

=  0.00000001 

JSR  =  60  dB 

Notch  Width  =  2 

JSR  =  45  dB 

^  =  0.000000001 

Threshold  =13.0 

JSR  =  65  dB 

Notch  Width  =  2 

^  =  0.000000001 

JSR  =  50  dB 

JSR  =  70  dB 

Threshold  =  17.0 

^1  =  0.0000000001 

Notch  Width  =  2 

JSR  =  75dB 

JSR  =  55  dB 

=  0.0000000001 

Threshold  =  20.0 

JSR  =  80  dB 

Notch  Width  =  2 

=  0.00000000001 

JSR  =  60  dB 

Threshold  =  18.0 

Notch  Width  =  4 

JSR  =  65dB 

Threshold  =  11.0 

Notch  Width  =  4 

JSR  =  70dB 

Threshold  =  15.0 

Notch  Width  =  4 

JSR  =  75dB 

Threshold  =15.0 

Notch  Width  =  4 

JSR  =  80  dB 

Threshold  =  25.0 

Notch  Width  =4 

Table  4.2:  Jammer  1-  Corr  /  MLT-TDED  -  Parameters  Used  for  Each  Value  of  JSR 
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In  figures  4.5  and  4.6,  the  curves  representing  BER  vs.  jammer  frequency  are  shown,  for  the 
different  receiver  structures  and  for  the  raw  error,  where  no  interference  suppression  technique  is 
incorporated.  The  plots  only  show  a  sampling  of  jammer  frequencies  with  points  taken  at  /  =  0.08  Hz, 
/  =  0.1//z,  /  =  0.12//z,  /  =  0.15//z,  /  =  0.18//z,  f  =  02Hz,  f  =  0234  Hz, 
f  =  025  Hz,  f  =  0.3  Hz,  f  =  035  Hz,  f  =  0.4  Hz,  and  /  =  0.45  Hz .  If  more  frequencies 
were  tested,  the  results  would  show  mwe  structure.  For  these  simulations,  the  default  parameters  used  are 
a  value  of  E/,  I  Ng  =6  dB  fw  the  channel,  and  a  value  of  JSR  =  20  dB  for  the  jammer.  Magnitudes 
of  /X  =  0.0000001  fw  the  predictive  and  /i  =  0.00001  for  the  two-sided  are  used,  while  a 
threshold  =  15.0  and  notch  width  =  1  are  employed  for  the  MLT-transform  domain  exciser.  The 
convergence  parameter  used  fw  the  adaptive  correlator  is/t  =  0.000001,  and  threshold  =  25.0  and 
notch  width  =  1  for  the  MLT-TDED  structure.  It  can  be  observed  that  only  the  MLT-TDED  system  is 
slightly  dependent  on  jammer  frequency. 


Figure  4.5 :  BER  vs.  Jammer  Frequency  for  Jammer  1  -  Filters  /  MLT-TDE 


JAMMER  1  -  Fixed  Tone 


Figure  4.6  :  BER  vs.  Jammer  Frequency  for  Jammer  1  •  Corr  /  MLT-TDED 


46 


4.2-  Results  for  Jammer  2 
4.2.1-  Swept  Tone 

The  interference  used  in  the  following  section  is  the  swept  tone  jammer,  where  a  triangular  wave 
is  used  to  vary  the  tone  frequency. 

In  Figure  4.7  and  4.8,  BER  vs.  Eh/No  for  the  different  receiver  systems  can  be  observed.  Also,  the 
performance  of  BPSK  baseband  signal  in  AWGN,  and  the  results  related  with  the  raw  error  are 
incorpwated  in  the  plots.  For  these  simulations,  the  default  parameters  used  for  the  swept  jammer  are 
JSR  =  20  dB  ,  a  normalized  frequency  /  =  0.25  Hz  and  a  sweep  frequency  =  0.00001  Hz . 

A  value  of  fJ.  =  0.00001  fw  the  predictive  and  fJ,  =  0.0001  for  the  two-sided  filter  are  used,  while  the 
MLT-transform  domain  excision  process  is  performed  with  threshold  =  13.0  and  notch  width  =  1. 
Values  of  11  =  0.00001  for  the  correlator,  and  threshold  =  15.0  and  notch  width  =  1  for  the  MLT- 

TDED  system  are  selected.  All  these  values  are  optimized  in  order  to  obtain  the  best  results  for  the 
systems.  It  can  be  noticed  from  Figure  4.7  that  the  adaptive  LMS  filters  and  the  MLT-TDE  structure 
introduce  improvements,  removing  the  interference  and  decreasing  the  BER.  The  predictive  filter  seems  to 
perform  slightly  worse  than  the  others.  Figure  4.8  shows  that  the  performance  of  the  MLT-TDED  system 
is  pretty  similar  to  the  ones  obtained  with  the  two-sided  filter  and  MLT-TDE  receiver.  On  another  hand, 
the  adaptive  correlatOT  presents  poor  results  in  comparison  to  the  other  systems. 


Figure  4.7 :  BER  vs.  Eb/No  for  Jammer  2-Triang  -  Filters  /  MLT-TDE 
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Figure  4.8  ;  BER  vs.  Eb/No  for  Jammer  2-Triang  -  Corr  /  MLT-TDED 


Figures  4.9  and  4.10  present  the  plots  of  BER  vs.  JSR  for  the  different  structures,  where  a  value 
of  Ef,  I N I,  =  6  dB  is  used  for  the  channel.  The  default  parameters  used  for  the  jammer  are  normalized 
frequency  of  /  =  0J25  Hz  and  sweep  frequency  of  f„„p  =  0.00001  Hz .  Tables  4.3  and  4.4  show  the 
optimized  values  for  the  convergence  factcv  fl  of  the  adaptive  filters  and  correlator,  and  exciser 
parameters  of  the  MLT  systems;  the  optimization  is  done  for  each  JSR  value.  From  the  graphs,  it  can  be 
concluded  that  none  of  the  receiver  systems  is  able  to  efficiently  operate  when  exposed  to  this  jammer. 
After  values  of  JSR  ^20  — 35  dB ,  the  perfOTmance  of  the  various  structures  starts  presenting  poor 
results,  even  tough  the  MLT-TDED  system  seems  to  handle  the  interference  slightly  better. 
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Figure  4.10  :  BER  vs.  JSR  for  Jammer  2-Triang  -  Corr  /  MLT-TDED 
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Table  4.3:  Jammer  2  Triangle  -  Filters  /  MLT-TDE  -  Parameters  Used  for  Each  Value  of  JSR 
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ADAPTIVE  CORRELATOR 


JSR=  10  dB 
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JSR=  15  dB 
^1  =  0.0001 
JSR  =  20dB 
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JSR  =  25dB 
^  =  0.00001 
JSR  =  30dB 

=  0.000001 
JSR=35dB 
^  =  0.000001 
JSR  =  40dB 
4  =  0.000001 
JSR  =  45dB 
4  =  0.0000001 
JSR  =  50dB 
4  =  0.0000001 
JSR  =  55dB 

4  =  0.00000001 
JSR  =  60dB 

4  =  0.00000001 
JSR  =  65dB 
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4  =  0.0000000001 
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Threshold  =  11.0 
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JSR  =  20  dB 

Threshold  =  15.0 
Notch  Width  =  1 
JSR  =  25dB 

Threshold  =11.0 
Notch  Width  =  2 
JSR=30dB 

Threshold  =  11.0 
Notch  Width  =  2 
JSR  =  35  dB 

Threshold  =11.0 
Notch  Width  =  2 
JSR  =  40  dB 

Threshold  =  11.0 
Notch  Width  =  2 
JSR  =  45  dB 

Threshold  =  11.0 
Notch  Width  =  4 
JSR  =  50dB 

Threshold  =  13.0 
Notch  Width  =  4 
JSR  =  55  dB 

Threshold  =11.0 
Notch  Width  =  4 
JSR=60dB 

Threshold  =  14.0 
Notch  Width  =  4 
JSR  =  65  dB 

Threshold  =  18.0 
Notch  Width  =  6 
JSR  =  70  dB 

Threshold  =  22.0 
Notch  Width  =  6 


JSR  =  75  dB 

Threshold  =  22.0 
Notch  Width  =  6 
JSR  =80  dB 

Threshold  =  30.0 

_ _  Notch  Width  =  6 

Table  4.4:  Jammer  2  Triangle  -  Corr  /  MLT-TDED  -  Parameters  Used  for  Each  Value  of  JSR 
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Figures  4.11  to  4.15  present  the  results  obtained  for  BER  vs.  sweep  frequency  for  the  different 
receiver  structures,  using  several  values  fOT  the  parameters  of  the  interference  suppression  techniques, 
namely  /X  ot  the  threshold  and  notch  width .  For  these  simulations,  the  default  parameters  used  are 
JSR  =  20  dB  ,  a  normalized  frequency  of  /  =  0.25  Hz ,  and  Ef,  I =  6  dB .  The  MLT  systems 
present  good  results  and  are  relatively  insensitive  to  sweep  frequency.  It  can  also  be  observed  that  their 
performance  is  rather  independent  of  the  parameters  used  fix'  the  excision  process.  On  another  hand,  the 
LMS  filters  and  correlator  have  hard  time  dealing  with  the  jamming  effects  after  values  of 

/sweep  ~  Hz .  Apparently,  the  LMS  algorithm  is  unable  to  track  the  jammer  when  its  frequency 

changes  too  rapidly.  Besides,  the  performance  of  these  systems  relies  heavily  on  the  convergence 
parameter  jj,  utilized;  the  higher  the  sweep  frequency,  the  lower  the  parameter  needs  to  be  in  order  to 
reduce  the  mean-square  enor  and  remove  the  jamming. 


Figure  4.11 :  BER  vs.  Sweep  Frequency  for  Jammer  2-Triang  •  1-S  Filter 
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Figure  4.12  :  BER  vs.  Sweep  Frequency  for  Jammer  l-Triang  -  2-S  Filter 


Figure  4.13  :  BER  vs.  Sweep  Frequency  for  Jammer  2-Triang  -  MLT-TDE 
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Figure  4.14 :  BER  vs.  Sweep  Frequency  fw  Jammer  2-Triang  -  Corr 


Figure  4.15 :  BER  vs.  Sweep  Frequency  for  Jammer  2-Triang  -  MLT-TDED 
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4.2.2-  Chirp  Jamming 

In  this  section,  the  chirp  jamming  is  added  to  the  channel  in  order  to  produce  the  interference. 
Here,  a  ramp  wave  is  used  to  vary  the  tone  frequency. 

Figures  4.16  and  4.17  show  the  BER  vs.  Eb/No  for  the  four  receiver  structures,  as  well  as  for  the 
results  obtained  with  BPSK  signal  in  AWGN  and  raw  errcff.  The  default  parameters  used  for  the  jammer 
are  JSR  =  20  dB ,  a  normalized  frequency  /  =  0.25  Hz  and  a  sweep  frequency 
/swap  —  0.00001  Hz .  The  optimized  values  for  the  system  parameters  are  ^  =  0.00001  for  the 
predictive  and  jU  =  0.0001  Cw  the  two-sided  filtCT,  threshold  =  1 1.0  and  notch  width  =  1  for  the 
MLT-TDE  system,  a  value  of  /X  =  0.00001  for  the  correlator,  and  threshold  =  10.0  and 
notch  width  =  1  for  the  MLT-TDED  system.  It  can  be  noticed  from  the  figures  that  all  the  techniques 
perform  efficiently  in  removing  the  jammer,  and  their  results  are  pretty  similar,  except  for  the  predictive 
filter  whose  performance  is  slightly  inferior  than  the  others. 


Figure  4.16 :  BER  vs.  Eb/No  for  Jammer  2-Ramp  -  Filters  /  MLT-TDE 


Figure  4.17 :  BER  vs.  Eb/No  for  Jammer  2-Ramp  -  Corr  /  MLT-TDED 


Figures  4.18  and  4.19  show  BER  vs.  JSR  fw  the  different  receiver  systems,  using  a  value  of 
!  N  ^  —  6dB  fcH’  the  Gaussian  noise  in  the  channel.  The  default  parameters  used  for  the  jammer  are 
ncrmalized  frequency  of  /  =  025  Hz  and  sweep  frequency  of  =  0.00001  Hz  .  Tables  4.5  and 

4.6  show  the  tqttimized  values  fw  the  systems  parameters  selected  fw  each  value  of  JSR.  As  the  plots 
indicate,  the  four  receiver  structures  show  pow  performance  in  removing  the  chirp  jamming.  After  values 
of  JSR  ~  30  dB,  the  results  start  deteriorating  rapidly,  making  an  exception  for  of  the  MLT-TDED 
system,  whose  performance  is  slightly  better. 
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Figure  4.19  :  BER  vs.  JSR  for  Jammer  2-Ramp  -  Corr  /  MLT-TDED 
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Threshold  =  43.0 
Notch  Width  =  3 
JSR  =  60dB 

Threshold  =  48.0 
Notch  Width  =  3 
JSR  =65  dB 

Threshold  =  50.0 
Notch  Width  =  0 
JSR  =  70dB 

Threshold  =  47.0 
Notch  Width  =  0 
JSR  =  75  dB 

Threshold  =  45.0 
Notch  Width  =  1 
JSR  =  80  dB 

Threshold  =  80.0 
Notch  Width  =  2 


Table  4.5:  Jammer  2  Ramp  -  Filters  /  MLT-TDE  -  Parameters  Used  for  Each  Value  of  JSR 


ADAPTIVE  CORRELATOR 


MLT-TDED 


JSR=  10  dB 
4  =  0.00001 
JSR=  15  dB 
4  =  0.00001 
JSR  =  20dB 
4  =  0.00001 
JSR  =  25dB 
4  =  0.000001 
JSR  =  30dB 
4  =  0.000001 
JSR  =  35  dB 
4  =  0.000001 
JSR  =  40dB 
4  =  0.0000001 
JSR  =  45dB 
4  =  0.0000001 
JSR  =  50dB 

4  =  0.00000001 
JSR  =  55  dB 

4  =  0.00000001 
JSR  =  60dB 

4  =  0.00000001 
JSR  =  65dB 

4  =  0.000000001 
JSR  =  70  dB 

4  =  0.000000001 
JSR  =  75  dB 

4  =  0.0000000001 
JSR=80dB 

4  =  0.0000000001 


JSR  =  10  dB 

Threshold  =  10.0 
Notch  Width  =  1 
JSR=  15  dB 

Threshold  =  10.0 
Notch  Width  =  1 
JSR  =  20dB 

Threshold=  10.0 
Notch  Width  =  1 
JSR  =  25dB 

Threshold  =  14.0 
Notch  Width  =  2 
JSR  =  30dB 

Threshold  =  12.0 
Notch  Width  =  2 
JSR  =  35  dB 

Threshold  =15.0 
Notch  Width  =  2 
JSR  =  40  dB 

Threshold  =11.0 
Notch  Width  =  2 
JSR  =  45  dB 

Threshold  =11.0 
Notch  Width  =  2 
JSR  =  50dB 

Threshold  =  10.0 
Notch  Width  =  2 
JSR=55dB 

Threshold  =13.0 
Notch  Width  =  4 
JSR  =  60  dB 

Threshold  =13.0 
Notch  Width  =  4 
JSR  =  65  dB 

Threshold  =  17.0 
Notch  Width  =  4 
JSR  =  70dB 

Threshold  =  20.0 
Notch  Width  =  4 
JSR  =  75dB 

Threshold  =  20.0 
Notch  Width  =  4 
JSR  =80  dB 

Threshold  =  24.0 
Notch  Width  =  4 


Table  4.6:  Jammer  2  Ramp  -  Corr  /  MLT-TDED  -  Parameters  Used  for  Each  Value  of  JSR 
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Figures  4.20  to  4.24  show  BER  vs.  sweep  frequency  fw  the  various  receiver  structures,  using 
many  different  values  fw  //  and  the  exciser  parameters.  The  default  parameters  used  are  JSR  —  20dB 
and  a  normalized  frequency  of  /  =  0.25  Hz  for  the  jammer,  and  a  value  of  E,,  I =  6dB  for  the 
channel.  The  MLT  systems  provide  a  significant  improvement  over  the  no-suppression  case,  along  with 
no  real  dependence  on  the  sweep  frequency  and  excisitm  process  parameters.  On  the  contrary,  the  LMS 

algorithm  is  unable  to  track  the  jammer  after  values  of  ftwtq>  “10  *  Hz  ,  and  its  performance  is  closely 
tied  to  the  value  of  the  convergence  parameter  fl . 


Figure  4.20 ;  BER  vs.  Sweep  Frequency  for  Jammer  2-Ramp  -  1-S  Filter 


60 


Figure  4.21 :  BER  vs.  Sweep  Frequency  for  Jammer  2-Ramp  -  2-S  Filter 


Figure  4.22  :  BER  vs.  Sweep  Frequency  for  Jammer  2-Ramp  -  MLT-TDE 


Figure  4.23  :  BER  vs.  Sweep  Frequency  for  Jammer  2-Ramp  -  Corr 


Figure  4.24  :  BER  vs.  Sweep  Frequency  for  Jammer  2-Ramp  -  MLT-TDED 


JAMMER  2  RAMP  -  Chirp  Jamming 


10’  -  Raw  Error:  -- 

•  Thres=1 0.0  /  Notch  Wictth=1 : 
Thres=8.0  /  Notch  Wicith=1 :  ... 

;  Thres=10.0  /  Notch  Width=2:  _ 

Thres=1 0.0  /  Notch  Width=0: 


Code  Order=6  /  Eb/No=6.0  /  JSR=20.0  /  Jammer  Freq=0.25 
Vector  Length:  MLT-TDED=64 


Sweep  Frequency 


4.3-  Results  for  Jammer  3 

This  section  shows  the  results  obtained  for  each  of  the  communication  receiver  systems  in 
presence  of  Gaussian  noise  jamming. 

Figure  4.25  shows  BER  vs.  Eh/N,,  fw  the  filters  and  MLT-based  transform  domain  exciser,  along 
with  the  results  for  BPSK  signal  in  AWGN  and  raw  error.  The  jammer  parameters  used  are 
JSR  =  20  dB  ,  a  normalized  frequency  /  =  0.234  Hz  ,  and  a  jammer  bandwidth  of  BW  =  0.1  Hz . 
A  value  of  /i  =  0.000001  fw  the  predictive  and  /X  =  0.0001  for  the  two-sided  filter  are  used,  while  a 
threshold  =  9.0  and  notch  width  =  2  are  employed  fcr  the  MLT-transform  domain  exciser.  The 
two-sided  filter  and  MLT-TDE  system  present  good  and  similar  interference  suppression  capabilities;  in 
contrast,  fo  higher  values  of  I *  3  tiS ,  the  predictive  filter  presents  poor  performance,  much 
infericff  than  the  one  obtained  with  the  two  previous  structures. 
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Figure  4.26  shows  the  plots  of  BER  vs.  Eh/No  fOT  the  LMS  correlator  and  MLT-TDED  system. 
Fot  these  simulatitxis,  the  default  parameters  used  for  the  jammer  are  JSR  =  20  dB  ,  a  normalized 
frequency  f  =  0.234  Hz ,  and  a  jammer  bandwidth  of  BW  =  0.1  Hz .  A  value  of  /I  =  0.000001  for 
the  correlator,  and  threshold  =  9.0  and  notch  width  =  1  for  the  MLT-TDED  system  are  used.  Note 
that  the  COTrelatw  structure  presents  similar  results  to  the  two-sided  filter  and  MLT-TDE  system  as  shown 
in  Figure  4.25,  while  the  perfOTmance  of  the  MLT-TDED  structure  is  slightly  inferior,  particularly  after 
the  value  of  I ~SdB. 


JAMMER  3  -  NB  Gaussian 


Figure  4.26 :  BER  vs,  Eb/No  for  Jammer  3  -  Corr  /  MLT-TDED 


In  Figures  4.27  and  4.28,  the  plots  of  BER  vs.  JSR  fw  the  different  structures  can  be  seen.  The 
values  used  fw  the  jammer  are  =  6  c£S  ,  a  nOTmalized  frequency  of  /  =  0.234  Hz  and  a 

jammer  bandwidth  of  BW  =  0.1  Hz .  The  optimized  parameters  of  the  interference  suppression 
techniques  fw  each  value  of  JSR  are  shown  in  Tables  4.7  and  4.8.  After  values  of  JSR  ~  30  dB ,  all  the 
different  structures  show  pow  performance  in  dealing  with  the  interference;  although  the  two-sided  LMS 
filter  and  the  MLT-TDED  structure  provide  a  slightly  better  performance  than  the  other  receivers. 
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Figure  4.27 :  BER  vs.  JSR  for  Jammer  3  -  Filters  /  MLT-TDE 


Figure  4.28  :  BER  vs.  JSR  for  Jammer  3  -  Corr  /  MLT-TDED 
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PREDICTIVE  FILTER 
JSR=10dB 

H  =  0.000001 

JSR  =  15  dB 
^  =  0.0000001 
JSR  =  20dB 
4  =  0.000001 
JSR  =  25dB 
4  =  0.000001 
JSR  =  30  dB 
4  =  0.0000001 
JSR  =  35dB 
4  =  0.000001 
JSR  =  40dB 
4  =  0.0000001 
JSR  =  45dB 
4  =  0.0000001 
JSR  =  50dB 
4  =  0.0000001 
JSR=55dB 
4  =  0.0000001 
JSR  =  60  dB 
4  =  0.0000001 
JSR  =  65dB 

4  =  0.0000001 
JSR  =  70dB 
4  =  0.0000001 
JSR  =  75dB 
4  =  0.0000001 
JSR  =  80dB 

4  =  0.00000001 


TWO-SIDED  FILTER 

JSR=  10  dB 
4  =  0.0001 
JSR  =  15  dB 
4  =  0.00001 
JSR  =  20dB 
4  =  0.0001 
JSR  =  25dB 
4  =  0.000001 
JSR  =  30dB 
4  =  0.000001 
JSR=35dB 
4=0.000001 
JSR  =  40dB 
4  =  0.0000001 
JSR  =  45dB 
4  =  0.0000001 
JSR  =  50dB 
4  =  0.0000001 
JSR  =  55dB 

4  =  0.00000001 
JSR  =  60dB 

4  =  0.00000001 
JSR  =  65dB 

4  =  0.00000001 
JSR  =  70dB 

4  =  0.000000001 
JSR  =  75  dB 

4  =  0.000000001 
JSR  =  80dB 

4  =  0.0000000001 


JSR  = 


JSR  = 


MLT-TDE 
:  10  dB 

Threshold  =  9.0 
Notch  Width  =  2 
:  15  dB 

Threshold  =  9.0 
Notch  Width  =  2 
20  dB 

Threshold  =  9.0 
Notch  Width  =  2 
25  dB 

Threshold  =11.0 
Notch  Width  =  2 
30  dB 

Threshold  =11.0 
Notch  Width  =  2 
35  dB 

Threshold  =  15.0 
Notch  Width  =  2 
40  dB 

Threshold  =  19.0 
Notch  Width  =  2 
45  dB 

Threshold  =  20.0 
Notch  Width  =  3 
50  dB 

Threshold  =  24.0 
Notch  Width  =  3 
55  dB 

Threshold  =  35.0 
Notch  Width  =  3 
60  dB 

Threshold  =  45.0 
Notch  Width  =  4 
65  dB 

Threshold  =  47.0 
Notch  Width  =  4 
70  dB 

Threshold  =  47.0 
Notch  Width  =  5 
75  dB 

Threshold  =  47.0 
Notch  Width  =  5 
80  dB 

Threshold  =  49.0 
Notch  Width  =  5 
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ADAPTIVE  CORRELATOR 

MLT-TDED 

JSR=  10  dB 

JSR=  10  dB 

11  =  0.00001 

Threshold  =  9.0 

JSR=  15  dB 

Notch  Width  =  4 

|i  =  0.000001 

JSR=15dB 

JSR  =  20  dB 

Threshold  =  9.0 

H  =  0.000001 

Notch  Width  =  4 

JSR  =  25  dB 

JSR  =  20dB 

11  =  0.000001 

Threshold  =  9.0 

JSR=30dB 

Notch  Width  =  1 

11  =  0.000001 

JSR  =  25dB 

JSR  =  35dB 

Threshold  =  9.0 

[i  =  0.000001 

Notch  Width  =  4 

JSR  =  40dB 

JSR=30dB 

\i  =  0.000001 

Threshold  =  9.0 

JSR  =  45dB 

Notch  Width  =  4 

H  =  0.0000001 

JSR  =  35dB 

JSR  =  50dB 

11  =  0.0000001 

Threshold  =  9.0 
Notch  Width  =  4 

JSR  =  55dB 

JSK.  —  40  dii 

Threshold  =  12.0 

11  =  0.0000001 

Notch  Width  =  4 

JSR  =  60dB 

JSR  =  45dB 

11  =  0.0000001 

Threshold  =  13.0 

JSR  =  65  dB 

Notch  Width  =  4 

11  =  0.0000001 

JSR  =  50dB 

JSR  =  70dB 

Threshold  =  15.0 

11  =  0.00000001 

Notch  Width  =  4 

JSR  =  75dB 

JSR  =  55dB 

11  =  0.00000001 

Threshold  =  15.0 

JSR=80dB 

Notch  Width  =  4 

11  =  0.00000001 

JSR  =  60dB 

Threshold  =  20.0 
Notch  Width  =  5 

JSR  =  65  dB 

Threshold  =  28.0 
Notch  Width  =  5 

JSR  =  70dB 

Threshold  =  35.0 
Notch  Width  =  5 

JSR  =  75dB 

Threshold  =  40.0 
Notch  Width  =  6 

JSR  =  80  dB 

Threshold  =  40.0 

Notch  Width  =  6 
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Figures  4.29  and  4.30  show  the  curves  representing  BER  vs.  jammer  frequency  Pot  the  different 
receiver  structures  and  for  the  raw  error.  The  plots  only  show  a  sampling  of  jammer  frequencies  with 
points  taken  at  /  =  0.08  Hz,  /  =  01  Hz,  f  =  012  Hz,  f  =  015  Hz,  f  =  018  Hz, 
f  =  02  Hz,  f  =  0234  Hz,  f  =  025  Hz,  f  =  0.3  Hz,  f  =  035  Hz,  f  =  0.4  Hz,  and 
/  =  0.45  Hz .  For  these  simulations,  the  default  parameters  used  for  the  jammer  are 
Ef,  I =  6  dB  ,  JSR  =  20  dB  and  a  jammer  bandwidth  of  BW  =  0.1  Hz .  A  value  of 
^  =  0.000001  fw  the  predictive  and  ^  =  0.0001  for  the  two-sided  filter  are  used,  while  a 
threshold  =  9.0  and  notch  width  =  2  are  employed  for  the  MLT-transfotm  domain  exciser.  Also,  a 
value  of  11  =  0.000001  fw  the  correlatw,  and  threshold  =  9.0  and  notch  width  =  1  for  the  MLT- 
TDED  system  are  utilized  to  perform  the  simulations.  The  results  show  that  the  LMS  filtCTS  and  correlator 
are  somewhat  dependent  on  jammer  frequency. 


JAMMER  3  •  NB  Gaussian 


Jammer  Frequency 


Figure  4.29  :  BER  vs.  Jammer  Frequency  fw  Jammer  3  -  Filters  /  MLT-TDE 
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JAMMER  3  -  NB  Gaussian 


Figure  4.30  :  BER  vs.  Jammer  Frequency  for  Jammer  3  -  Corr  /  MLT-TDED 


In  Figures  4.31  and  4.32,  the  plots  of  BER  vs.  jammer  bandwidth  can  be  seen  for  the  different 
receiver  systems  and  for  the  raw  error  case.  The  curves  only  show  a  sampling  of  jammer  bandwidths  with 
results  obtained  at  BW  =  0.03  Hz,  BW  =  0.05  Hz,  BW  =  01  Hz,  BW  =  0.15  Hz, 
BW=02Hz,  BW=015Hz,  BW  =  03  Hz,  BW  =  035  Hz,  BW  =  OA  Hz,  and 
BW  =  0.45  Hz.  The  default  parameters  are  !  N  „=6dB  ,  JSR  =  20  dB  and  a  normalized 

frequency  of  /  =  0.234  Hz .  The  parameters  fOT  the  interference  suppression  techniques  include  a  value 
of  /i  =  0.000001  for  the  predictive  and  fX  =  0.0001  for  the  two-sided  filter,  a  threshold  =  9.0  and 
notch  width  =  2  for  the  MLT-transform  domain  exciser,  a  value  of  fl  =  0.000001  for  the  correlator, 
and  values  of  threshold  =  9.0  and  notch  width  =  1  for  the  MLT-TDED  system.  For  all  the 
structures,  the  results  deteriorate  as  the  jammer  bandwidth  increases:  the  predictive  filter  shows  inferior 
performance,  while  the  two-sided  filter  and  the  MLT-TDED  system  seem  to  tolerate  larger  bandwidth 
slightly  better. 
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Figure  4J1 :  BER  vs.  Jammer  Bandwidth  for  Jammer  3  •  Filters  /  MLT-TDE 


Figure  4.32  ;  BER  vs.  Jammer  Bandwidth  for  Jammer  3  -  Corr  /  MLT-TDED 
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4.4-  Results  for  Jammer  4 
4.4.1-  Pulsed  Tone 

This  section  presents  the  results  achieved  for  each  of  the  receivers  when  the  channel  faces  pulsed 
tone  jamming. 

In  figure  4.33,  the  plot  of  HER  vs.  Eh/No  for  the  one  and  two-sided  LMS  filters,  and  for  the  MLT- 
based  transform  domain  exciser  is  shown.  Besides,  the  results  obtained  with  BPSK  baseband  signal  in 
AWGN  alone,  along  with  the  performance  obtained  without  any  interference  suppression  technique,  are 
presented  in  the  same  curve.  For  these  simulations,  the  default  parameters  used  for  the  jammer  are 
JSR  =  20  dB  and  a  normalized  frequency  /  =  0.1  Hz ,  while  values  of  juk  —  0.000001  Hz 

and  Duty  Cycle =  50  %  are  utilized  fw  the  masking  function.  A  magnitude  of  /X  =  0.00001 
for  the  predictive  and  /X  =  0.0001  for  the  two-sided  filter  are  used,  while  a  threshold  =  12.0  and 
notch  width  =  1  are  utilized  for  the  MLT-TDE  structure.  As  can  be  seen  in  the  figure,  the  three 
suppression  structures  introduce  significant  improvement  over  the  no-suppression  case,  although  the 
predictive  filter  presents  a  slighdy  inferior  performance  than  the  other  two  systems. 
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Figure  4.34  shows  BER  vs.  Eb/No  fw  the  adaptive  LMS  correlator  and  MLT-TDED  system.  Also, 
the  performance  of  BPSK  baseband  signal  in  AWGN  alone,  as  well  as  the  performance  of  the  raw  error 
are  added.  Fca-  these  simulations,  the  default  parameters  used  for  the  jammer  are  JSR  =  20  dB  and  a 
normalized  frequency  f=  0.1  Hz,  while  values  of  =0.000001  Hz  and 

Duty  =  50  %  are  utilized  few  the  masking  function,  jl  =  0.000001  for  the  correlator, 

and  threshold  =  15.0  and  notch  width  =  2  few  the  MLT-TDED  system  are  used.  These  parameters 
were  e^Wiroized  in  ewder  to  achieve  best  perfewmance  few  the  receivers.  Both  structures  present  similar 
results,  which  are  closer  to  the  AWGN-ewily  performance  than  the  two-sided  filter  and  MLT-TDE 
structures  presented  in  Figure  4.33. 


Figure  4.34  :  BER  vs.  Eb/No  for  Jammer  4-Tone  -  Corr  /  MLT-TDED 
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Figures  4.35  and  4.36  show  the  curves  of  HER  vs.  JSR  for  the  different  structures.  The 
convergence  factor  fX  fw  the  adaptive  filters  and  correlator,  and  exciser  parameters  for  the  MLT  systems 
are  optimized  for  each  value  of  JSR,  and  can  be  seen  in  Tables  4.9  and  4.10.  A  value  of 
Eu  I N  =6  dB  is  used  for  the  Gaussian  noise  included  in  the  channel.  The  default  parameters  used  for 

the  jammer  are  a  normalized  frequency  /  =  0.1  Hz  ,  and  values  of  ^  =  0.000001  Hz  and 

Duty  Cycle^ii  =  50  %  fw  the  masking  function.  It  can  be  seen  that  the  LMS  filters  provide  the 

best  perfcffmance,  in  particular  the  two-sided  filter,  showing  moderate  dependence  of  BER  on  JSR  .  All 
the  other  systems,  provide  poor  performance  after  values  of  JSR  ~  35  dB ,  even  though  the  MLT-TDED 

system  seems  to  reasonably  reduce  the  jamming  effects  up  to  JSR  ~  65  dB . 
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JAMMER  4  -  Tone  w/  Mask  Func 


Figure  4.36  :  BER  vs.  JSR  for  Jammer  4-Tone  -  Corr  /  MLT-TDED 
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PREDICTIVE  FILTER 

JSR=  10  dB 
^1  =  0.00001 
JSR=  15  dB 
^1  =  0.00001 
JSR  =  20  dB 
4  =  0.00001 
JSR  =  25dB 
4  =  0.00001 
JSR  =  30dB 
4  =  0.000001 
JSR  =  35dB 
4  =  0.000001 
JSR  =  40dB 
4  =  0.0000001 
JSR  =  45  dB 
4  =  0.0000001 
JSR=50dB 

4  =  0.0000001 
JSR  =  55dB 
4  =  0.0000001 
JSR  =  60  dB 

4  =  0.00000001 
JSR  =  65dB 

4  =  0.00000001 
JSR  =  70dB 

4  =  0.000000001 
JSR  =  75  dB 

4  =  0.000000001 
JSR  =  80dB 

4  =  0.0000000001 


two-sided  filter 

JSR=10dB 
4  =  0.0001 
JSR=  15  dB 
4  =  0.0001 
JSR  =  20  dB 
4  =  0.0001 
JSR  =  25dB 
4  =  0.0001 
JSR  =  30  dB 
4  =  0.00001 
JSR  =  35dB 
4  =  0.00001 
JSR  =  40dB 
4  =  0.000001 
JSR  =  45dB 
4  =  0.000001 
JSR  =  50dB 
4  =  0.0000001 
JSR  =  55  dB 
4  =  0.0000001 
JSR  =  60dB 

4  =  0.00000001 
JSR  =  65dB 

4  =  0.00000001 
JSR  =  70dB 

4  =  0.000000001 
JSR  =  75dB 

4  =  0.000000001 
JSR  =  80dB 

4  =  0.0000000001 


MLT-TDE 


JSR=  10  dB 

Threshold  =  9.0 
Notch  Width  =  1 
JSR=  15  dB 

Threshold  =  12.0 
Notch  Width  =  1 
JSR  =  20dB 

Threshold  =  12.0 
Notch  Width  =  1 
JSR  =  25dB 

Threshold  =  12.0 
Notch  Width  =  1 
JSR  =  30dB 

Threshold  =  12.0 
Notch  Width  =  1 
JSR=35dB 

Threshold  =  16.0 
Notch  Width  =  2 
JSR  =  40dB 

Threshold  =  20.0 
Notch  Width  =  2 
JSR  =  45dB 

Threshold  =  18.0 
Notch  Width  =  4 
JSR  =  50dB 

Threshold  =  18.0 
Notch  Width  =  5 
JSR  =  55dB 

Threshold  =  18.0 
Notch  Width  =  5 
JSR  =  60  dB 

Threshold  =  22.0 
Notch  Width  =  6 
JSR  =  65  dB 

Threshold  =  22.0 
Notch  Width  =  7 
JSR  =  70dB 

Threshold  =  24.0 
Notch  Width  =  7 
JSR  =  75dB 

Threshold  =  24.0 
Notch  Width  =  7 
JSR=80dB 

Threshold  =  24.0 
Notch  Width  =  7 


Table  4.9:  Jammer  4  Tone  -  Filters  /  MLT-TDE  -  Parameters  Used  for  Each  Value  of  JSR 
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ADAPTIVE  CORRELATOR 


MLT-TDED 


JSR  =  10  dB 

=  0.000001 
JSR  =  15  dB 
M.  =  0.000001 
JSR  =  20  dB 

\i  =  0.000001 

JSR  =  25  dB 
^  =  0.000001 
JSR  =  30  dB 

H  =  0.000001 

JSR=35dB 

ji  =  0.000001 

JSR  =  40dB 
4  =  0.000001 
JSR  =  45  dB 
4  =  0.0000001 
JSR  =  50dB 

4  =  0.00000001 
JSR  =  55  dB 

4  =  0.000000001 
JSR  =  60  dB 

4  =  0.0000000001 
JSR  =  65  dB 

4  =  0.0000000001 
JSR  =  70  dB 

4  =  0.0000000001 
JSR  =  75dB 

4  =  0.000000001 
JSR  =  80  dB 

4  =  0.000000001 


JSR  =  10  dB 

Threshold  =  11.0 
Notch  Width  =  1 
JSR=  15  dB 

Threshold  =  15.0 
Notch  Width  =  2 
JSR  =  20  dB 

Threshold  =  15.0 
Notch  Width  =  2 
JSR  =  25  dB 

Threshold  =  15.0 
Notch  Width  =  2 
JSR  =  30dB 

Threshold  =  17.0 
Notch  Width  =  2 
JSR  =  35dB 

Threshold  =  21.0 
Notch  Width  =  2 
JSR  =  40dB 

Threshold  =  23.0 
Notch  Width  =  3 
JSR  =  45dB 

Threshold  =  25.0 
Notch  Width  =  4 
JSR  =  50  dB 

Threshold  =  27.0 
Notch  Width  =  4 
JSR  =  55  dB 

Threshold  =  27.0 
Notch  Width  =  4 
JSR  =  60dB 

Threshold  =  27.0 
Notch  Width  =  5 
JSR  =  65  dB 

Threshold  =  35.0 
Notch  Width  =  5 
JSR  =  70dB 

Threshold  =  35.0 
Notch  Width  =  6 
JSR  =  75  dB 

Threshold  =  35.0 
Notch  Width  =  6 
JSR  =  80  dB 

Threshold  =  37.0 
Notch  Width  =  6 


Table  4.10:  Jammer  4  Tmie  -  Corr  /  MLT-TDED  -  Parameters  Used  for  Each  Value  of  JSR 
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Figures  4.37  and  4.38  present  BER  vs.  jammer  frequency,  for  the  different  receiver  structures  and 
for  the  raw  error.  The  curves  only  show  a  sampling  of  jammer  frequencies  with  points  taken  at 
/  =  0.08//z,  f  =  QXHz,  f  =  Q12Hz.  /  =  0.15//z,  f  =  0X%Hz,  f  =  02Hz, 
f  =  0234  Hz,  f  =  025  Hz,  f  =  0.3  Hz,  f  =  035  Hz,  f  =  0.4  Hz,  and  /  =  0.45  Hz . 
For  these  simulations,  the  default  parameters  used  are  values  of  I =  6(IB  fw  the  channel, 
JSR  =  20  dB  fw  the  jammo-,  and  /,„*  yj,*.  =  0.000001  Hz  and  Duty  Cycle ^  =50%  for 
the  masking  function.  The  parameters  used  for  the  interference  suppression  techniques  are  =  0.00001 
for  the  predictive  and  /X  =  0.0001  fw  the  two-sided  filter,  threshold  =  12.0  and  notch  width  =  1 
for  the  MLT-TDE  structure,  fl  =  0.000001  fcff  the  correlatOT,  and  threshold  =  15.0  and 
notch  width  =  2  fw  the  MLT-TDED  system.  It  can  be  observed  that  all  the  structures  are  able  to 
efficiently  remove  the  jamming  effects,  showing  no  real  dependence  on  jammer  frequency. 


JAMMER  4  -  Tone  w/  Mask  Func 


Figure  4.37 :  BER  vs.  Jammer  Frequency  for  Jammer  4-Tone  -  Filters  /  MLT-TDE 
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JAMMER  4  -  Tone  w/  Mask  Func 


Jammer  Frequency 


Figure  438 :  BER  vs.  Jammer  Frequency  fw  Jammer  4-Tone  -  Corr  /  MLT-TDED 


4.4.2-  Pulsed  Gaussian  Jamming 

This  section  shows  the  results  obtained  fm*  each  of  the  communication  receiver  systems  in 
presence  of  pulsed  Gaussian  noise  jamming  in  the  channel. 

In  Figure  4.39,  the  plot  of  BER  vs.  Eh/No  for  the  filters  and  MLT-based  transform  domain 
exciser,  as  well  as  the  results  for  BPSK  signal  in  AWGN  and  raw  error  are  presented.  The  jammer 
parametCTS  used  are  75/?  =  20  dB  ,  a  normalized  frequency  of  /  =  0.1  Hz ,  and  a  jammer  bandwidth 
of  BW  =  0.05  Hz,  while  the  values  for  the  masking  function  are  =  0.000001 //z  and 

Duty  =  50  %  .  A  value  of  (I  =  0.00001  for  the  predictive  and  =  0.00001  for  the 

two-sided  filter  are  used,  while  a  threshold  =  10.0  and  notch  width  =  1  are  employed  for  the  MLT- 
TDE  structure.  The  two-sided  filter  and  MLT-TDE  system  present  reaswiable  and  somewhat  similar 
interference  suppression  capabilities;  on  another  hand,  for  higher  values  of  =  4  d!S ,  the 

predictive  filter  presents  poOT  perfOTmance,  with  higher  values  for  the  BER  than  the  other  two  structures. 
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Figure  4.39  :  BER  vs.  Eb/No  for  Jammer  4-NB  -  Filters  /  MLT-TDE 


Figure  4.40  plots  BER  vs.  Eb/No  for  the  LMS  correlatw  and  MLT-TDED  system.  For  these 
simulations,  the  default  parameters  used  fw  the  jammer  are  JSR  =  20  dB  ,  a  normalized  frequency 
/  =  0.1  Hz  ,  and  a  Jammer  bandwidth  of  BW  =  0.05  Hz ,  while  the  default  values  utilized  to  represent 
the  masking  function  are  =0.000001  Hz  and  Duty  Cycle =50%.  A  value  of 

fj.  =  0.00001  fw  the  correlator,  and  threshold  =  13.0  and  notch  width  =  2  for  the  MLT-TDED 

system  are  used.  Both  structures  produce  results  that  are  similar  to  those  for  the  two-sided  filter  and  MLT- 
TDE  system,  even  though  the  results  of  the  adaptive  correlator  slightly  deteriorate  after 
EJN,^5dB. 
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Figures  4.41  and  4.42  show  the  curves  representing  BER  vs.  JSR  for  the  different  structures. 
Values  of  I  =  6dB  ,  UOTmalized  frequency  of  /  =  0.1  Hz  ,  and  jammer  bandwidth  of 
BW  =  0.05  Hz  are  used  fw  the  jammer,  while  for  the  masking  function  =  0.000001  Hz 

and  Duty  Cycle^,^  ;s«c  =  50  %  are  employed.  The  optimized  parameters  of  the  interference 
suppression  techniques  fw  each  value  of  JSR  are  shown  in  Tables  4. 1 1  and  4. 12.  As  can  be  observed  from 
the  curves,  after  values  of  JSR  =  30  — 40  dB  all  the  different  structures  show  pow  performance  in 

improving  the  BER;  although  the  two-sided  LMS  filter  and  the  MLT-TDED  structure  present  better 
results  than  the  other  receivers. 
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Figure  4.41 :  BER  vs.  JSR  for  Jammer  4-NB  •  Filters  /  MLT-TDE 


Figure  4.42  :  BER  vs.  JSR  for  Jammer  4-NB  -  Corr  /  MLT-TDED 
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PREDICTIVE  FILTER 

TWO-SIDED  FILTER 

JSR  =  10  dB 

JSR  =  10  dB 

ti  =  0.00001 

11  =  0.00001 

JSR  =  15  dB 

JSR  =  15  dB 

|i  =  0.00001 

11  =  0.00001 

JSR  =  20  dB 

JSR  =  20dB 

11  =  0.00001 

11  =  0.00001 

JSR  =  25dB 

JSR  =  25dB 

11  =  0.00001 

11  =  0.00001 

JSR=30dB 

JSR  =  30  dB 

11  =  0.00001 

11  =  0.00001 

JSR  =  35dB 

JSR  =  35dB 

11  =  0.000001 

11  =  0.000001 

JSR  =  40dB 

JSR  =  40  dB 

11  =  0.000001 

11  =  0.000001 

JSR  =  45dB 

JSR  =  45dB 

11  =  0.0000001 

11  =  0.0000001 

JSR  =  50dB 

JSR  =  50dB 

11  =  0.0000001 

It  =  0.0000001 

JSR  =  55  dB 

JSR  =  55  dB 

11  =  0.00000001 

11  =  0.00000001 

JSR  =  60  dB 

JSR  =  60dB 

11  =  0.00000001 

It  =  0.00000001 

JSR  =  65dB 

JSR  =  65dB 

11  =  0.000000001 

11  =  0.000000001 

JSR  =  70dB 

JSR  =  70dB 

11  =  0.0000000001 

11  =  0.000000001 

JSR  =  75dB 

JSR  =  75dB 

11  =  0.000000000001 

11  =  0.0000000001 

JSR=80dB 

JSR=80dB 

11  =  0.000000000001 

11  =  0.0000000001 

MLT-TDE 

JSR  =  10  dB 

Threshold  =  10.0 
Notch  Width  =  1 
JSR=  15  dB 

Threshold  =11.0 
Notch  Width  =  2 
JSR  =  20dB 

Threshold  =  10.0 
Notch  Width  =  1 
JSR  =  25dB 

Threshold  =11.0 
Notch  Width  =  2 
JSR  =  30  dB 

Threshold  =11.0 
Notch  Width  =  2 
JSR  =  35dB 

Threshold  =  13.0 
Notch  Width  =  2 
JSR  =  40  dB 

Threshold  =15.0 
Notch  Width  =  2 
JSR  =  45  dB 

Threshold  =  17.0 
Notch  Width  =  4 
JSR  =  50  dB 

Threshold  =  21.0 
Notch  Width  =  4 
JSR  =  55dB 

Threshold  =  21.0 
Notch  Width  =  5 
JSR  =  60  dB 

Threshold  =  21.0 
Notch  Width  =  5 
JSR  =  65dB 

Threshold  =  25.0 
Notch  Width  =  6 
JSR  =  70  dB 

Threshold  =  29.0 
Notch  Width  =  7 
JSR  =  75dB 

Threshold  =  29.0 
Notch  Width  =  8 
JSR  =  80dB 

Threshold  =  29.0 
Notch  Width  =  8 
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ADAPTIVE  CORRELATOR 

JSR=  10  dB 

p.  =  0.00001 

JSR=  15  dB 

H  =  0.00001 

JSR  =  20dB 

H  =  0.00001 

JSR  =  25dB 
^  =  0.000001 
JSR=30dB 
^  =  0.000001 
JSR  =  35dB 
4  =  0.000001 
JSR  =  40dB 
4  =  0.000001 
JSR  =  45dB 
4  =  0.0000001 
JSR=50dB 
4  =  0.0000001 
JSR=55dB 
4  =  0.0000001 
JSR  =  60dB 
4  =  0.0000001 
JSR  =  65dB 

4  =  0.00000001 
JSR  =  70dB 

4  =  0.00000001 
JSR=75dB 

4  =  0.00000001 
JSR  =  80dB 

4  =  0.00000001 


Table  4.12:  Jammer  4  Noise  -  Corr  /  MLT-TDED  -  Parameters  Used  for  Each  Value  of  JSR 


MLT-TDED 

JSR=  10  dB 

Threshold  =  9.0 
Notch  Width  =  1 
JSR=  15  dB 

Threshold  =  13.0 
Notch  Width  =  2 
JSR  =  20dB 

Threshold  =  13.0 
Notch  Width  =  2 
JSR  =  25dB 

Threshold  =  14.0 
Notch  Width  =  2 
JSR=30dB 

Threshold  =  14.0 
Notch  Width  =  2 
JSR  =  35dB 

Threshold  =  14.0 
Notch  Width  =  2 
JSR  =  40dB 

Threshold  =  14.0 
Notch  Width  =  2 
JSR  =  45  dB 

Threshold  =  14.0 
Notch  Width  =  2 
JSR=50dB 

Threshold  =  14.0 
Notch  Width  =  2 
JSR=55dB 

Threshold  =  14.0 
Notch  Width  =  2 
JSR  =  60dB 

Threshold  =  14.0 
Notch  Width  =  5 
JSR  =  65dB 

Threshold  =  14.0 
Notch  Width  =  5 
JSR  =  70dB 

Threshold  =  14.0 
Notch  Width  =  6 
JSR  =  75  dB 

Threshold  =  14.0 
Notch  Width  =  6 
JSR  =  80  dB 

Threshold  =  14.0 
Notch  Width  =  6 


83 


Figures  4.43  and  4.44  show  the  curves  representing  BER  vs.  jammer  frequency  for  the  different 
receiver  structures  and  fw  the  raw  errw.  The  plots  only  show  a  sampling  of  jammer  frequencies  with 
points  taken  at  /  =  0.08  Hz,  /  =  OJ  Hz,  f  =  0J2  Hz,  f  =  0J5  Hz,  /  =  018  Hz, 
f  =  02Hz,  f  =  0234  Hz,  f  =  025  Hz,  f  =  0.3  Hz,  f  =  035  Hz,  f  =  0.4  Hz,  and 
/  =  0.45  Hz .  Fot  these  simulations,  the  default  parameters  used  for  the  jammer  are 
!  N g  =  6dB  ,  JSR  =  20 dB  and  a  jammer  bandwidth  of  BW  =  0.05  Hz ,  while  values  of 
fmsk  June  ~  0.000001  Hz  and  Duty  Cycle ys,*.  =  50  %  are  utilized  for  the  masking  frequency.  A 
value  oi  11  =  0.00001  for  the  predictive  and  fl  =  0.00001  for  the  two-sided  filter  are  used,  while  a 
threshold  =  10.0  and  notch  width  =  1  are  employed  for  the  MLT-TDE  structure.  Also,  a  value  of 
11  =  0.00001  fcM-  the  correlator,  and  threshold  =  13.0  and  notch  width  =  2  for  the  MLT-TDED 
system  are  utilized  to  perfcHm  the  simulations.  The  results  show  that  the  structures  are  ratho-  independent 
on  jamma*  frequency. 


Figure  4.43  :  BER  vs.  Jammer  Frequency  for  Jammer  4-NB  -  Filters  /  MLT-TDE 
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Figure  4.44 :  BER  vs.  Jammer  Frequency  for  Jammer  4-NB  -  Corr  /  MLT-TDED 


In  Figures  4.45  and  4.46.  the  plots  of  BER  vs.  jammer  bandwidth  can  be  seen  for  the  different 
receiver  systems  and  for  the  raw  error  case.  The  curves  only  show  a  sampling  of  jammer  bandwidths  with 
results  obtained  at  BW  =  0.03  Hz,  BW  =  0.05  Hz.  BW  =0.1  Hz, 

BW  =  015  Hz, 

BW=02Hz,  BW=025Hz,  BW=03Hz,  BW  =  0.35  Hz,  BW=0.AHz,  and 
BW  =  0.45  Hz .  The  default  parameters  fcff  the  jammer  otq  E,,  I =  6  dB  ,  JSR  =  20  dB  and  a 
normalized  frequency  of  /  =  0.1  Hz ;  ~  0.000001  Hz  and  Duty  Cycle fu„c  —50% 

are  used  fw  the  masking  frequency.  The  parameters  for  the  interference  suppression  techniques  include  a 
value  of  jU  =  0.00001  fw  the  predictive  and  ^  =  0.00001  for  the  two-sided  filter,  a 

threshold  =  10.0  and  notch  width  =  1  fw  the  MLT-TDE  structure,  a  value  of  jU  =  0.00001  for  the 
correlatOT,  and  threshold  =  13.0  and  notch  width  =  2  for  the  MLT-TDED  system.  In  all  the 
structures,  the  BER  increases  with  the  jammer  bandwidth;  the  predictive  filter  shows  inferior 
performance,  while  the  two-sided  filter  and  the  MLT-TDED  system  seem  to  somewhat  perform  better 
under  higher  jammer  bandwidth  scenarios. 
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Figure  4.45  ;  BER  vs.  Jammer  Bandwidth  for  Jammer  4-NB  -  Filters  /  MLT-TDE 


Figure  4.46 :  BER  vs.  Jammer  Bandwidth  for  Jammer  4-NB  -  Corr  /  MLT-TDED 
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CHAPTER  5 

Conclusions  and  Discussions 


5.1-  Conclusions 

The  performance,  in  terms  of  low  bit  error  rate,  of  the  different  receiver  systems  depends  heavily 
on  the  jamming  signal  and  parameters  used,  and  therefcffe  there  is  no  “ideal”  structure  in  the  sense  of  a 
particular  system  that  would  be  able  to  handle  any  type  of  situation.  However,  the  satisfactory  achievement 
of  the  LMS  filters  fw  the  fixed  tone,  and  of  the  MLT-TDED  system  for  the  sweep  and  Gaussian  noise 
Jammers  can  be  emphasized. 

5.1.1-  Adaptive  Filters 

The  simulation  results  indicate  that  LMS  adaptive  filtering  improves  bit  error  rate  perfwmance 
in  a  variety  of  narrowband  jammer  environments.  Nevertheless,  for  higher  values  of  JSR  ,  the  value  of 
the  H  parameter  has  to  be  lowered  down,  to  decrease  the  additional  distortion  introduced  by  the  jammer. 
The  adaptive  LMS  filter  perfoms  pretty  well  with  the  fixed  tone,  but  has  problems  keeping  track  of  the 
changes  involved  in  the  swept-tone  jammers.  Indeed,  since  the  jamm«  fi'equency  is  now  changing 
constantly  with  these  jammCTs,  the  uansfo'  function  of  the  adaptive  filter  must  also  be  changing 
continually  to  track  them.  Hence,  the  LMS  algwthm  must  not  only  converge,  but  must  converge  quickly 
enough  to  adjust  any  change  in  the  jammer;  thus,  the  convergence  factor  H  of  the  adaptive  filter  becomes 
an  important  parameter,  as  it  governs  the  rate  at  which  the  filter  converges.  As  mentioned  earlier  in  this 
report,  the  larger  the  value  of  /i  ,  the  faster  the  convergence.  Therefore,  it  is  necessary  to  maintain  a  large 
value  of  ^  ,  as  the  jammer  frequency  is  changing  quickly,  however,  a  larger  fj.  introduces  mwe  code 
distortion,  making  difficult  the  overall  performance  of  the  filter.  Moreover,  it  can  be  observed  that  the 
performance  relies  heavily  on  the  value  of  ^  used;  consequently,  careful  tuning  is  necessary  in  wder  for 
the  filters  to  q)erate  adequately. 

In  general,  the  two-sided  filter  performs  better  than  the  predictive,  and  requires  lower 
implementation  complexity  with  only  25  OpS  (Operation  per  Sample),  while  the  predictive  filter  needs 
28  OpS  .  It  can  therefore  be  concluded  that  the  two-sided  structure  is  much  more  efficient  than  the 
predictive. 

5.1.2-  MLT-Transform  Domain  Exciser 

The  MLT-TDE  system  performs  pretty  well  as  long  as  the  interference  power  remains  lower  than 
a  certain  limit  Essentially,  if  the  jammer-to-signal  ratio  is  high  enough,  the  windowing  effects  of  the 
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block  processing  may  cause  the  jammer  sidelobes  to  also  transcend  the  threshold,  and  as  a  result,  a  large 
portion  of  the  signal  energy  is  removed,  increasing  the  bit  error  rate. 

The  results  also  point  to  the  feet  that  the  MLT  systems  are  relatively  insensitive  to  sweep 
frequency,  and  that  their  perfOTmance  is  rather  independent  of  the  parameters  used  fw  the  excision 
process,  which  means  that  the  system  does  not  require  fine  tuning. 

5.1.3-  Adaptive  Correlator 

Regarding  the  tracking  performance  of  this  receiver,  the  adaptive  correlator  shows  the  highest 
convergence  time,  as  its  tracking  capabilities  are  unable  to  remove  the  interference  when  the  jammer 
sweeps  too  rapidly,  this  is  why  a  higher  number  of  symbols  to  wait  is  used  for  this  structure.  Indeed,  this 
receiver  involves  a  training  period  where  the  known  sequence  is  generated  in  the  receiver.  Moreover,  this 
structure  presents  poor  results  with  large  values  of  for  the  swept-tone  jammers. 

Ifowever,  the  adaptive  cwrelatOT  involves  much  less  computation  in  comparison  with  all  the 
other  receivers,  with  only  4  OpS  (Operation  per  Sample),  due  to  its  relatively  simple  structure. 

5.1.4-  MLT-Transform  Domain  Excision  and  Detection 

The  MLT-TDED  system  perfwms  pretty  well  as  long  as  the  jammer-to-signal  ratio  is  not  high 
enough  to  originate  signal  energy  destruction,  even  though  its  performance  is  pretty  satisfactory  compared 
with  the  othCT  receiver  structures,  specially  fw  the  swept-tone  and  Gaussian  jammers. 

Besides,  as  the  results  indicate,  the  tracking  capabilities  of  the  MLT-TDED  system  are  excellent, 
and  it  can  definitely  be  concluded  that  this  structure  presents  the  shortest  convergence  time. 

However,  regarding  the  implementation  complexity,  the  MLT-TDED  system  requires  29  OpS  , 
and  therefore  is  the  computationally  most  expensive  system.  Anyhow,  the  lapped  transfearm  domain 
excision  and  detection  algcxithm  must  be  considered  as  an  efficient  method  for  narrowband  interference 
suppressicxi  applications. 

5.1.5-  Additional  Comments 

It  can  be  noticed  from  Tables  4.1  to  4.12  that  the  higher  the  value  of  JSR  ,  the  lower  the  value 
of  p.  in  the  LMS  algorithm  to  compensate  fex  the  excess  power  of  the  jammer,  and  the  higher  the  values 

of  threshold  and/o-  notch  width  in  the  excision  process,  in  order  to  be  able  to  remove  larger  amcamts 
of  interference  energy.  Also,  the  values  of  p  used  for  the  predictive  filter  are  usually  equal  ot  lower  than 
the  ones  used  in  the  two-sided  structure.  Indeed,  the  results  obtained  with  the  first  structme  are  frequently 
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infericH'  than  the  ones  achieved  with  the  two-sided,  and  therefwe  a  lower  values  of  /i  has  to  be  used  in 
ordCT  to  reduce  the  excess  mean-square  enor. 

It  seems  like  the  overall  perfainance  of  the  various  receiver  systems  presents  better  outcomes 
with  the  fixed  tone  than  with  the  pulsed  tone;  apparently,  the  changing  shape  of  the  jammer  produces 
higher  error  rate.  HbwevCT,  better  performance  is  obtained  with  the  pulsed  Gaussian  noise  than  with  the 
uninterrupted  version  of  it.  This  result  might  be  due  to  the  fact  that  a  value  of  BW  =  0.05  Hz  is  used 
fra:  the  pulsed  noise,  while  a  higher  value  of  BW  =  0.1  Hz  is  applied  to  the  non-pulsed  noise  jammer. 
Also,  the  results  obtained  with  chirp  jamming  are  slightly  better  than  the  ones  with  swept  tone.  This 
might  be  due  to  the  fwm  of  the  triangular  wave  itself,  which  exhibits  more  alterations  than  the  ramp 
wave,  and  therefore  is  harder  to  track. 

5.2-  Future  Research 

Some  simulations  were  run  for  the  filters  and  correlate  receivers  using  the  RLS  algorithm  for 
the  adaptive  filths.  The  RLS  algcxithm  is  much  mwe  computationally  expensive  than  the  LMS,  even 
though  it  is  supposed  to  converge  faster.  However,  as  fw  the  different  systems  studied,  the  performance  of 
the  RLS  algCKithm  was  somewhat  similar  to  the  LMS  algorithm’s.  Also,  nonlinear  LMS  and  RLS  filters 
could  be  analyzed  [12]. 

LMS  and  RLS  adaptive  filtering  in  the  frequency  domain  can  give  significant  reductions  in 
computation  over  the  conventional  time-domain  approach  [4],  and  could  also  be  incorporated  as  another 
communication  receiver  structure.  However,  this  system  has  the  drawback  that  it  may  not  be  able  to  react 
as  quickly  to  changes  in  the  input  as  the  simple  exciser. 

Additional  processing  could  be  incorporated  in  the  different  systems,  like  compensation  for  the 
adaptive  filter,  and  errOT-correcting  codes  combined  with  interleaving. 
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